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ABSTRACT 


Data  are  given  in  the  form  of  attenuation  factors  for  the  expo¬ 
sure  due  to  gamma  radiation  transmitted  by  a  ribbed  slab.  The  ribbed  slab 
is  made  of  concrete  and  is  similar  to  one  which  has  been  used  in  experi¬ 
mental  studies  conducted  at  the  University  of  Illinois.  The  source  radi¬ 
ation  was  assumed  to  be  that  of  Co-60  with  source  spectrum  degradation  due 
to  the  self-shielding  of  the  source.  Four  angles  of  incidence,  0°,  45°, 
60°,  and  75°,  were  considered.  In  addition,  the  effect  of  a  beam  of  radi¬ 
ation  incident  with  directions  diverging  2.5°  on  either  side  of  45°  was 
studied  in  a  rather  crude  fashion.  Attenuation  factors  for  1.25  MeV  gamma 
radiation  incident  normally  on  a  simulated  wood  floor  are  included  in  an 
appendix. 
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I.  INTRODUCTION 

The  primary  purpose  of  the  calculations  described  in  this  report 
is  to  provide  a  set  of  theoretical  data  to  be  compared  with  results  ob¬ 
tained  in  ribbed-slab  shielding  experiments  which  have  been  conducted 
at  the  University  of  Illinois.  Preliminary  comparison  between  experi¬ 
mental  results  and  the  data  of  this  report  indicates-  that  satisfactory 
agreement  can  be  expected.  However,  the  analysis  of  the  experimental 
data  is  not  yet  complete  and  detailed  comparisons  will  be  made  in  a 
later  report. 

The  basic  geometry  which  was  considered  is  a  plane,  horizontal,  homo¬ 
geneous  slab,  with  ribs  constructed  from  the  same  material  resting  par¬ 
allel  to  each  other  on  the  slab.  A  plane,  monodirectional  source  of 
gamma  radiation  is  assumed  to  be  incident  on  the  side  of  the  slab 
opposite  the  ribs.  Detectors  are  placed  at  various  distances  from  the 
ribbed  side  of  the  slab  and  at  various  horizontal  positions  relative 
to  the  ribs.  The  responses  of  these  detectors  are  expressed  as  atten¬ 
uation  factors.  In  this  report,  the  attenuation  factor  for  a  given 
detector  is  defined  as  the  ratio  of  the  exposure  received  by  the  de¬ 
tector  when  the  ribbed  slab  is  present  to  the  exposure  it  would  re¬ 
ceive  when  the  ribbed  slab  is  absent. 

The  second  section  of  this  report  describes  the  Monte  Carlo  cal¬ 
culation  of  the  contribution  of  scattered  radiation  to  the  attenuation 
factor  and  an  analytical  calculation  of  the  contribi'tion  of  uncoil'  ed 
radiation.  In  the  third  section,  data  comparable  to  the  experimental 
data  mentioned  above  are  presented  and  discussed.  Attenuation  factors 
for  radiation  incident  normally  on  a  simulated  wood  floor  are  tabu¬ 
lated  in  Appendix  A.  The  co:  .->uter  program  which  was  used  for  the  cal¬ 
culations  is  described  in  Appendix  B. 

II.  CALCULATION  OF  ATTENUATION  FACTORS 

A.  General  Description  of  the  Problem 

Figure  1  shows  a  cross-section  view  of  a  ribbed  slab.  The  ribs 
are  infinitely  long  and  perpend* cular  to  the  plane  of  the  figure.  The 
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z-axis  is  perpendicular  to  the  ribs.  The  y-axis  (not  shown  in  Fig¬ 
ure  1)  is  parallel  to  the  ribs.  Photon  directions  are  specified  by 

direction  cosines,  u  ,  u  ,  and  u  ,  relative  to  the  coordinate  axes. 

x  y  z 

The  initial  photon  direction  is  defined  by  the  cosine  of  the  polar 

angle  0  and  the  azimuthal  angle  <D  so  that  the  initial  direction 
o  o 

cosines  are  given  by 


u  = 

sin0 

COSCD  , 

xo 

o 

o 

u  = 

sin0 

sincp  , 

yo 

o 

o 

u  = 
zo 

COS0 

o 

• 

The  attenuation  factor  is  computed  as  a  function  of  detector 
height  above  the  slab  and  the  horizontal  position  of  the  detector 
relative  to  the  ribs.  It  is  designated  by  A  (0  ,H,X)  where, 


r 


Here,  (x,  z  )  are  the  coordinates  of  the  detector  position,  z  is 
d  s 

the  slab  thickness  not  including  the  ribs,  and  x^  is  half  the  sepa¬ 
ration  distance  between  rib  centers  (see  Figure  1).  Because  the 
rib  structure  extends  to  x  =  ±  °  ,  and  the  ribs  are  infinitely  long, 

A  (0  ,H,X)  is  a  periodic  function  of  X  .  Thus,  it  was  calculated 
f  o 

only  for  0  £  X  S  1  .  The  dependence  of  A  (0  ,HtX)  on  co  ,  the  slab 

i  o  o 

material,  and  the  rib  configurations  has  not  been  indicated  explicitly 
because  the  calculations  in  this  report  were  done  only  for  coq  =  0 
and  for  a  single  slab  material  and  slab-rib  geometry. 

The  attenuation  factors  A  (0  ,H,X)  were  computed  in  two  steps. 

f  o 

The  contribution  due  to  scattered  radiation  was  calculated  by  Monte 
Carlo  methods  and  the  contribution  of  unscattered  radiation  was 
determined  analytically.  Thus,  the  attenuation  factor  may  be  written 
as  the  sum  of  two  terms, 

A  (0  ,H,X)  =  A*(0  , H, X)  +  A°(0  H,X) 
i  o  I  o  i  o 
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where  the  superscript  s  refers  to  scattered  radiation  and  o 
refers  to  unscattered  radiation. 

B.  Monte  Carlo  Calculation  for  Scattered  Radiation 

Source  points  were  selected  by  dividing  the  interval  -x^  £  x  £  x 

into  a  certain  number  of  increments  and  choosing  the  midpoint  of  each 

increment  as  the  starting  point.  An  identical  number  of  histories 

originated  at  each  point.  The  initial  photon  direction  was  the  same 

for  each  history.  The  source  energy  spectrum  was  assumed  to  consist 

of  N  discrete  energies  E  each  having  probability  w  .  Thus, 

N 

l  ”»  =  1  ' 

n  =  1 

At  the  beginning  of  each  history,  one  of  the  energies  E  was 

n 

selected  at  random  according  to  the  probability  function  P(E  )  =  w  . 

n  n 

Only  photon  trajectories  within  the  region  -x  S  x  £  x  (see 

r  r 

Figure  1)  were  considered.  If  the  photon  left  this  region  by  crossing 

the  boundary  x  =  x  (x  =  ~x  ),  it  was  treated  as  a  photon  entering 
r  r 

the  region  at  x  =  -x  (x  =  x  ).  Each  time  the  position  for  an  inter- 

r  r 

action  was  determined,  the  coordinates  (x,z)  of  the  interaction  point 

were  checked  to  see  if  -x  <  x  <  x  and  z  <  z  <  z  .  If  this 

s  s  -  s  r 

was  found  to  be  the  case,  the  photon  was  moved  along  its  trajectory 

until  it  re-entered  the  ribbed  slab  or  until  z  became  equal  to  z  . 

r 

If  the  photon  re-entered  the  slab,  a  position  for  a  new  interaction 
point  was  selected. 

The  Monte  Carlo  calculation  was  designed  to  accommodate  any 
photon  energy  below  10  MeV.  Thus,  the  photoelectric  effect,  Compton 
scattering,  and  pair  production  were  taken  into  account.  At  the  be¬ 
ginning  of  each  photon  history,  the  photon  was  assigned  a  statistical 
weight  of  unity.  Each  time  the  photon  had  a  collision,  this  statis¬ 
tical  weight  was  multiplied  by  the  probability  that  the  interaction 
was  not  a  photoelectric  absorption.  Then  either  a  Compton-scattering 
interaction  or  a  pair- production  interaction  was  selected  using  the 
appropriate  probabilities  for  these  interactions,  given  that  photo- 
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electric  absorption  did  not  occur.  A.  new  photon  energy  and  direction 
were  determined  and  the  history  continued.  If  a  pair-production  inter¬ 
action  was  selected,  the  statistical  weight  was  doubled,  the  new  photon 
energy  was  set  equal  to  0.511  MeV,  and  the  new  photon  direction  was 
sampled  from  an  isotropic  distribution.  If  a  Compton-scattering  inter¬ 
action  was  selected,  the  new  photon  energy  and  direction  were  sampled 
from  the  Klein-Nishina  distribution  for  unpolarized  photons.  A  photon 
history  was  terminated  when  the  photon  was  transmitted  or  reflected 
by  the  slab,  or  when  its  energy  dropped  below  a  specified  cutoff 
energy. 

The  interval  -x  £  x  £  x  was  divided  into  a  number  of  detec- 
r  r 

tion  intervals  as  illustrated  in  Figure  1.  When  a  photon  was  trans¬ 
mitted  by  the  slab,  it  crossed  the  detector  plane  z  =  z,  by  pass- 
ing  through  one  of  these  detection  intervals.  Scores  for  several 
detector  planes  were  recorded  simultaneously.  The  contribution  to 
flux  was  estimated  by  dividing  the  statistical  weight  of  the  photon 

by  the  direction  cosine  u  of  the  photon  trajectory.  To  avoid  the 

z 

problem  of  the  infinite  variance  of  this  flux  estimator,  a  cutoff 

value  for  u  was  introduced  as  suggested  by  Clark1.  For  the  data 
z 

in  this  report,  contributions  to  the  flux  were  not  recorded  when 
the  direction  cosine  u^  of  the  trajectory  of  the  transmitted  pho¬ 
ton  was  less  than  0.01.  The  error  introduced  by  the  use  of  this 
cutoff  was  much  smaller  than  the  statistical  standard  deviation  of 
the  final  results.  The  average  exposure  was  calculated  for  each 
detection  interval  and  recorded  as  a  function  of  the  midpoint  of 
the  interval. 

C.  Analytical  Calculation  for  Uncollided  Radiation 

Attenuation  factors  for  the  uncollided  radiation  were  computed 
for  a  detector  located  at  the  midpoint  of  each  detection  interval 
used  in  the  Monte  Carlo  calculation.  This  part  of  the  calculation 
was  done  analytically.  The  main  complication  arose  from  the  fact 
that  the  photon  path  length  within  the  ribbed  slab  depended  on  both 
the  height  of  the  detector  above  the  slab  and  the  horizbntal  posi¬ 
tion  of  the  detector  with  respect  to  the  ribs.  Nevertheless,  the 
evaluation  of  this  path  length  was  straightforward  although  some¬ 
what  tedious  and  will  not  be  described  in  detail. 

Once  the  path  length  through  the  ribbed  slab  was  evaluated 
the  attenuation  factor  for  unc.ollided  radiation  was  calculated  using 
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the  formula, 

N 

y  W  H  (E  )  E  exp{-|i(E  )  t<0  H,X)} 

i-i  1  n  en  n  n  n  o 

Af(eo»H.x)  =  - , 

y  w  |1  (E  )  E 
Z_i  ,  n  en  n  n 

n  =  1 

where  p,gn  is  the  linear  energy  absorption  coefficient  for  air,  p,  is 
the  total  linear  attenuation  coefficient  for  the  ribbed  slab,  and 
t(0Q,H,X)  is  the  path  length  within  the  ribbed  slab. 

III.  ATTENUATION  FACTORS  FOR  A  PARTICULAR  RIBBED-SLAB  CONFIGURATION 
A.  Description  of  Input  Data 

In  this  section,  calculated  results  are  given  for  a  ribbed  slab 
whose  dimensions  correspond  closely  to  the  dimensions  of  the  ribbed 
slab  used  in  the  experimental  studies  at  the  University  of  Illinois. 
Referring  to  the  symbols  as  defined  in  Figure  1,  some  of  the  impor¬ 
tant  dimensions  were: 

x  =  6  1/8  inches 
r 

x  =  4  1/8  inches 
s 

z  =  4  inches 
s 

z  =  10  inches 
r 

The  slab  and  ribs  were  assumed  to  be  so-called  NBS  concrete 

3 

with  density  2.38  g^cm  .  The  composition  assumed  for  the  concrete 
is  listed  in  Table  1.  The  mass  interaction  coefficients  used  in 
the  calculation  are  listed  in  Table  2.  These  are  based  on  the 

2 

atomic  interaction  coefficient  data  used  by  Hubbell  and  Berger  . 

Also  listed  in  Table  2  are  the  mass  energy  absorption  coefficients 

for  air  which  were  used  in  the  calculation;  the  latter  are  also 

2 

from  Hubbell  and  Berger  . 

The  spectrum  of  gamma  radiation  emitted  by  the  experimental 

3 

source  was  calculated  in  an  earlier  report  .  When  the  source  was 

in  use,  a  1/8  inch  lead  filter  was  placed  over  the  end  of  the 

source  and  an  auxiliary  collimator  was  placed  above  the  filter. 

Hence, the  spectrum  for  the  calculation  was  taken  as  the  normalized 

spectrum  for  a  1/8  inch  filter  given  in  Table  A2  of  the  report 
3 

mentioned  above  .  This  spectrum  is  illustrated  in  rigure  2. 
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TABLE  1 


Composition  assumed  for  the  concrete  used  in  the 
ribbed  slafc^. 


Element 

Fraction 

by 

Weight 

H 

0.0056 

0 

0.4983 

Na 

0.0171 

Mg 

0.0024 

A1 

0.0456 

Si 

0.3158 

S 

0.0012 

K 

0.0192 

Ca 

0.0826 

Fe 

0.0122 

Sum 

1.0000 
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TABLE  2 

Mass  interaction  coefficients  for  concrete.  The  density  assumed 

3 

for  the  concrete  was  2. 38g/cm  .  Also  given  are  the  mass  energy 

absorption  coefficients  for  air.  Units  for  all  the  data  are 

2. 

cm  /g. 


£ 

Concrete 

Air 

Compton 

Pair 

Total 

Mass  Energy 

1  1 

Scattering 

Production 

Attenuation 

Absorption 

1  1 

Coefficient 

Coefficient 

Coefficient 

Coefficient 

0.010 

0.193 

0 

26.5 

4.61 

0.015 

0.190 

0 

8.01 

1.27 

0.020 

0.186 

0 

3.45 

0.511 

0.030 

0.180 

0 

1.12 

0.148 

0.040 

0.174 

0 

0.559 

0. 0668 

0.050 

0.169 

0 

0.361 

0. 0406 

0.060 

0.164 

0 

0.273 

0. 0305 

0.080 

0.156 

0 

0.200 

0. 0243 

0.100 

0,148 

0 

0.170 

0. 0234 

0.150 

0.134 

0 

0.140 

0.0250 

0.200 

0.122 

0 

0.125 

0.0268 

0.300 

0.106 

0 

0.107 

0.0287 

0.400 

0.0954 

0 

0.0957 

0.0295 

0.500 

0.0871 

0 

0. 0873 

0.0296 

0.600 

0.  0806 

0 

0.0807 

0.  0295 

0.800 

0. 0708 

0 

0. 0708 

0.0289 

1.000 

0. 0636 

0 

0. 0637 

0.0278 

1.500 

0.0517 

0.000155 

0.0518 

0.0254 

2.  000 

0.0441 

0.000629 

0.0447 

0. 0234 

3.000 

0. 0347 

0.00179 

0.0365 

0.0205 

4.000 

0.0290 

0. 00292 

0.0319 

0.0186 

5.  000 

0.0250 

0. 00394 

0.0290 

0.0174 

6.000 

0,0221 

0.00485 

0.0270 

0.0164 

8.000 

0.0181 

0. 00641 

0.0245 

0.0152 

10. 000 

0.0154 

0. 00771 

0.0231 

0.0145 
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B.  Results  for  a  Perfectly  Collimated  Plano  Source 

Attenuation  factors  are  given  in  Tables  3,  4,  5,  and  6  for  angles 

of  incidence  9q  =  0°,  45°,  60°,  and  75°.  The  initial  azimuthal  angle 

cp  was  set  equal  to  zero  in  all  cases.  In  the  tables,  data  are  listed 
°  s 

separately  for  A  (9  ,H,X)  and  A  (9  ,H,X)  .  The  percent  statistical 
f  o  f  o 

s 

standard  deviation  for  A. (9  ,H,X)  was  nearly  constant  as  a  function 

f  o 

of  H  and  X  for  a  given  value  of  8  ;  an  average  value  for  this 

o 

quantity  is  given  at  the  beginning  of  each  table.  The  data  for  each 

of  the  angles  9=0°  and  60°  are  based  on  100,000  histories.  The 
o 

data  for  each  of  the  angles  8q  =  45  and  75  are  based  on  10,000 
histories.  The  time  required  to  process  100,000  histories  on  the  IBM- 
7094  at  the  University  of  Illinois  is  about  thirty  minutes. 

s 

The  attenuation  factor  A. (9  ,H,X)  is  plotted  in  Figures  3  and  4 

f  o 

o  o 

for  angles  of  incidence  9  -  C  and  60  respectively.  For  both  angles 

o 

of  incidence,  the  results  show  a  strong  dependence  on  the  horizontal 
detector  position  X  when  H  =  0.49,  the  height  of  the  ribs.  However, 
the  data  very  rapidly  lose  their  dependence  on  X  as  the  detector 
height  increases;  for  H  >  1  ,  the  data  are  essentially  independent  of 

X  . 

The  result  of  including  the  uncollided  radiation  in  the  atten¬ 
uation  factors  for  these  angles  of  incidence  is  shown  in  Figures  5  and 
6.  As  can  be  seen,  in  this  case  a  dependence  on  X  persists  at  all 
values  of  H  .  The  amplitude  of  this  dependence  is,  of  course,  a 
function  of  the  angle  of  incidence.  For  the  60°  case,  the  main  effect 
of  changing  H  is  to  change  the  value  of  X  where  the  attenuation  fac¬ 
tor  A„(0  ,H,X)  has  its  maximum, 
f  o 

Since,  as  illustrated  in  Figure  3  and  4,  the  attenuation  factor 
for  scattered  radiation  rapidly  loses  its  dependence  on  X  as  H  in¬ 
creases,  it  is  of  interest  to  examine  the  relationship  between  the 
average  attenuation  factor, 

A®(9  )  =f  A®(8  ,H,X)  dX 

f  °  i  f  ° 

for  the  ribbed  slab  and  the  attenuation  factor  for  a  plane  slab  having 
the  same  average  mass  thickness.  Monte  Carlo  calculations  were 
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table  3 

ATTENUATION  FACTORS  FOR  0oa  o  DEGREES.  THE  FRACTIONAL 
STANDARD  DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  3  PERCENT 
FOR  THE  ATTENUATION  FACTORS  DUE  TO  SCATTERED  RADIATION  ONLY , 
DATA  IN  THE  TABLES  SHOULD  BE  MULTIPLIED  BY  0 . 1  TO  GET  THE 
CORRECT  ATTENUATION  FACTOR. 


X 

*49 

.98 

1.47 

1.96 

2.94 

3.92 

4.90 

5.88 

SCATTERED  RADIATION 

.025 

1  .03 

1  .75 

1  .88 

1,78 

1.91 

1  .81 

1.93 

1.86 

.075 

1.11 

1  .68 

1.79 

1.75 

1.78 

1  .88 

1  .81 

1.92 

.125 

1  .31 

1  .81 

1.81 

1.78 

1.81 

X  .89 

1  .72 

1  .85 

.175 

1  .76 

1 .75 

1  .89 

1  .88 

1.85 

1  .77 

1  .85 

1.84 

.225 

2.  10 

1 . 87 

1.79- 

1.80 

1.82 

1  .87 

1  .79 

1.81 

.275 

2.18 

1.77 

1  .90 

1 .86 

1 .83 

1  .79 

1  .81 

1.81 

.325 

2.06 

2.01 

1  .89 

1 .81 

1 .94 

1  .83 

1  .85 

1.77 

375 

2.16 

1  .89 

1.81 

1  .85 

1.80 

1  .81 

1  .83 

1.83 

.425 

2.32 

1  .84 

1  .77 

1.84 

1.78 

1  .88 

1  .88 

1.80 

.475 

2.25 

1  .88 

1  .88 

1.76 

1.86 

1  .90 

1  .83 

1.85 

.525 

2.31 

2.09 

1.95 

2.00 

1 .78 

1  .93 

1  .87 

1.80 

.575 

2.26 

1.89 

1  .81 

1.85 

1 .83 

1  .91 

1  ,83 

1.71 

.625 

2.19 

1  .83 

1  .81 

1.78 

1.88 

1  .90 

1  .79 

1.91 

.675 

2.17 

1.94 

1  .83 

1.91 

1 .93 

1  .88 

1  .83 

1.86 

.725 

2.23 

1  .92 

1  .85 

1.78 

1.85 

1  .76 

1  .87 

1 .75 

,775 

1  .93 

1.79 

1  .86 

1,85 

1.78 

1  .77 

1  .88 

1.82 

,825 

1  .84 

1 .81 

1.80 

1  .93 

1 .72 

1  .88 

1  .87 

1.92 

.875 

1  .31 

1.75 

1  .78 

1.72 

1.84 

1  .79 

1  ,87 

1.80 

.925 

1.07 

1 .68 

1.86 

1  .90 

1 .79 

1  .78 

1  .81 

1.84 

.975 

1.09 

1 . 76 

1.77 

1  .87 

1  .91 

1  .69 

1 .86 

1*94 

SCATTERED 

PLUS  UNSCATTERED  RADIATION 

.025 

1  .34 

2.06 

2.19 

2.08 

2,22 

2.1  1 

2*24 

2i  1  7 

.075 

1  ,42 

1.99 

2.09 

2.06 

2.08 

2.18 

2*  1  1 

2*22 

.125 

1  .62 

2.  12 

2.  1  l 

2,08 

2.12 

2.19 

2,03 

2*  16 

.175 

4.23 

4.21 

4.35 

4.34 

4.31 

4.23 

4,31 

4.30 

.225 

4.56 

4.33 

4.25 

4.26 

4.28 

4.33 

4.25 

4o  27 

.275 

4.64 

4.23 

4.36 

4.32 

4.29 

4.25 

4.27 

4.27 

.325 

4.52 

4.47 

4.35 

4.27 

4.40 

4  .29 

4.31 

4.23 

.375 

4.62 

4.35 

4.27 

4,31 

4.26 

4  .27 

4,29 

4.30 

.425 

4.78 

4.30 

4.23 

4.30 

4.24 

4.34 

4.34 

4.26 

.475 

4.71 

4.34 

4,34 

4,22 

4.32 

4.36 

4,29 

4.31 

.525 

4.77 

4.55 

4.41 

4.46 

4.24 

4.39 

4.33 

4 .26 

.575 

4.72 

4.35 

4.27 

4.31 

4.29 

4.37 

4.29 

4.  17 

.625 

4.65 

4.29 

4.27 

4.24 

4,34 

4.36 

4.25 

4.37 

.675 

4.63 

4.40 

4,29 

4.37 

4.39 

4.35 

4.29 

4.32 

.725 

4.69 

4.38 

4.31 

4.24 

4.31 

4.22 

4.33 

4.21 

.775 

4.39 

4.25 

4.32 

4.31 

4.24 

4.23 

4.34 

4.28 

.825 

4.30 

4.27 

4,26 

4.39 

4.18 

4.34 

4.33 

4.38 

,875 

1  .62 

2.06 

2.09 

2,02 

2. 15 

2.09 

2.17 

2.10 

.925 

1  .38 

1.98 

2.16 

2.20 

2,10 

2.09 

2.  1  1 

2.15 

.975 

1 .40 

2.07 

2.07 

2.  18 

2.22 

1  .99 

2.  1  1 

2.25 
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TABLE  4 

ATTENUATION  FACTORS  FOR  0O=  45  DEGREES.  THE  FRACTIONAL 
STANDARD  DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  13  PERCENT 
FOR  THE  ATTENUATION  FACTORS  DUE  TO  SCATTERED  RADIATION  ONLY. 
DATA  IN  THE  TABLES  SHOULD  BE  MULTIPLIED  BY  0.1  TO  GET  THE 
CORRECT  ATTENUATION  FACTOR. 


vd 

.49 

.98 

1  .47 

1.96 

2.94 

3.92 

4.90 

5.88 

SCATTERED  RADIATION 

.025 

.97 

1 .09 

1.78 

1.08 

1  .30 

1  .36 

.96 

1.11 

.075 

.77 

1  .24 

1  .28 

1.18 

1  .43 

1  .38 

1  .27 

1.38 

.125 

.61 

1.33 

1  .25 

1.30 

1  .06 

1  .09 

1  .21 

1  .07 

•  175 

.84 

1.24 

1  .03 

1.09 

.89 

.99 

1.18 

1.02 

.225 

.85 

1.02 

1  .04 

1.23 

1  .01 

1  .46 

1  .53 

1.19 

.275 

.91 

1  .39 

1.21 

1.27 

1  .  16 

1  .34 

1.26 

1.51 

.325 

1  .31 

.99 

1  .01 

1.66 

.92 

1  .02 

1  .  10 

1.08 

.375 

1 .17 

1.01 

1  .24 

1.18 

1  .08 

1 .06 

.94 

1.25 

.425 

1  .43 

1.43 

1  .03 

1.35 

1 .34 

1  .43 

1.10 

.80 

.475 

1  .41 

1  .  15 

1.24 

1.06 

1.21 

1  .40 

1.18 

1.45 

.525 

1.40 

1.  17 

1.28 

.95 

1.02 

1.30 

1.12 

1  .06 

.575 

1  .48 

1.14 

1  .  19 

1.04 

1 .33 

1  .32 

1  .55 

1  .36 

.625 

1  .66 

1.19 

1.23 

1.31 

1  *06 

1  .16 

1  .22 

1  .26 

.675 

1.20 

1.22 

1  .22 

1.30 

1.25 

1 .18 

1.25 

1.16 

.725 

1.46 

1.45 

1.12 

1.40 

1 .06 

1 .10 

1.23 

1.44 

.775 

1  .63 

1.11 

.97 

1.17 

1.34 

1  .06 

1.03 

1.16 

.825 

1  .53 

1.34 

1,  14  ‘ 

1.15 

1 .25 

1  .16 

1  .22 

1.28 

,875 

1  .22 

1.20 

1.30 

.97 

1.52 

1  .23 

1,44 

1.23 

.925 

1 .15 

1.14 

1  .20 

1,16 

1,45 

1 .16 

1.30 

1.09 

.975 

1.11 

1  .25 

1  .36 

1.26 

1  .40 

.90 

1.04 

1.21 

SCATTERED 

PLUS  UNSCATTERED  RADIATION 

.025 

1  .42 

1.77 

2.17 

1.85 

2.  17 

2.34 

2,06 

2,36 

.075 

1.  10 

2.  15 

1.57 

2.21 

2.60 

2.70 

2.65 

2.76 

•  125 

.86 

2.57 

1.47 

2,69 

2,44 

2.47 

2.59 

2.45 

.  175 

1  .03 

2.62 

1.23 

2.47 

2.27 

2.37 

2.57 

2.40 

.225 

1  .05 

2.40 

1.23 

2.61 

2.39 

2.84 

2.91 

2.57 

.275 

1 .10 

2.78 

1.41 

2.65 

2,54 

2.65 

2.43 

2.54 

.325 

1  .50 

2.37 

1  .22 

2.90 

2.02 

1  .99 

1.96 

1.84 

.375 

1  .43 

2,04 

1.53 

2.09 

1 .89 

1  .78 

1.58 

1.82 

.425 

1  .78 

2.  19 

1  .43 

2.02 

1 .94 

1  .96 

1  .57 

1  .22 

.475 

1  .88 

1.72 

1  .77 

1 .56 

1.65 

1  .79 

1  .53 

1.76 

.525 

2.04 

1  .59 

2.00 

1.32 

1 .35 

1  .59 

1.38 

1.29 

.575 

2.35 

1.45 

2.16 

1 ,31 

1 .57 

1  .53 

1  .74 

1.55 

.625 

2.82 

1.42 

2.55 

1.51 

1 .26 

1  .36 

1  .41 

1.45 

.675 

2.58 

1.4  1 

2.60 

1 .50 

1.45 

1  .37 

1  .44 

1.36 

.725 

2.84 

1.65 

2.51 

1 .59 

1.26 

1  .29 

1 .42 

1  .65 

.775 

3.01 

1  .30 

2.35 

1.37 

1 .54 

1  .28 

1.27 

1.44 

.825 

2.91 

1 .54 

2.46 

1.39 

1.51 

1  .46 

1.56 

1.65 

.875 

2.32 

1  .48 

2.28 

1.28 

1.87 

1  .63 

1.89 

1.73 

.925 

1  .96 

1.52 

1.92 

1  ,58 

1 .93 

1  .70 

1.90 

1.78 

.975 

1.71 

1.76 

1.90 

1.83 

2.04 

1  .62 

1  .86 

2.  13 
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table  s 

ATTENUATION  FACTORS  FOR  0O=I  60  DEGREES.  THE  FRACTIONAL 
STANDARD  DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  5  PERCENT  ' 
FOR  THE  ATTENUATION  FACTORS  DUE  TO  SCATTERED  RAD  I  AT  I  ON' ONLy « 
DATA  IN  THE  TABLES  SHOULD  BE  MULTIPLIED  BY  0*1  TO  GET  THE 
CORRECT  ATTENUATION  FACTOR.  ‘  1 


X 

.49 

S 

.98 

1.47 

1.96 

2.94 

3.92 

4.90 

5*83 

SCATTERED  RADIATION 

.025 

•  44 

.72 

.68 

.71 

.68 

•  66 

*  60 

*66 

.075 

.42 

.63 

.69 

.64 

.66 

.63 

♦  60 

*64 

.  125 

.38 

.60 

.61 

.63 

.64 

.65 

♦  61 

*  66 

.  175 

.42 

.67 

.63 

.  65 

.62 

.66 

♦  67 

*  64 

.225 

.53 

.70 

.61 

.67 

.65 

•  68 

♦  65 

*60 

.275 

.53 

.62 

.66 

.61 

.63 

.63 

*65 

*64 

.325 

.59 

.63 

.68 

.63 

.64 

.63 

.63 

*  66 

.375 

.61 

.64 

.58 

.61 

.62 

.64 

.64 

*62 

.425 

.73 

.63 

.63 

.61 

.60 

.60 

♦  60 

♦  61 

.475 

.69 

.  59 

.64 

.67 

.69 

.63 

♦  63 

•  63 

.525 

•  81 

.60 

•  61 

.68 

.62 

.67 

.65 

•  60 

.575 

.78 

.67 

.65 

.60 

.59 

.61 

.69 

•  68 

.625 

.78 

.65 

•  61 

.60 

.70 

.62 

.61 

*58 

.675 

.80 

.63 

.65 

.64 

.  62 

.63 

.65 

.70 

.725 

.82 

.60 

.59 

•  65 

.61 

.59 

.66 

.67 

,775 

.87 

.62 

•  66 

.68 

.65 

.70 

.65 

•  63 

,825 

.84 

.64 

.65 

.65 

.61 

.64 

♦  58 

♦  59 

.875 

.67 

.67 

.60 

.58 

.65 

.62 

♦  65 

*68 

.925 

.59 

.66 

.66 

.63 

.  66 

.65 

♦  69 

♦  62 

.975 

.50 

.64 

.71 

.65 

.66 

.67 

♦  69 

*69 

scattered 

PLUS  UNSCATTERED  RADIATION 

.025 

,68 

.85 

.81 

.83 

.88 

.84 

*72 

*79 

.075 

.61 

.75 

.81 

.77 

.91 

.77 

♦  72 

*80 

.  125 

.53 

.72 

.74 

.75 

.90 

♦  77 

♦  73 

*87 

.  175 

.55 

.80 

.75 

.78 

.  88 

.78 

*79 

*90 

.225 

.  56 

.82 

.74 

.80 

.91 

•  80 

♦  78 

*86 

.275 

•  65 

.75 

.78 

.77 

.  86 

.76 

♦  78 

*90 

.325 

.72 

.76 

.80 

.83 

.82 

♦  75 

♦  76 

*92 

.375 

.73 

.76 

.71 

•  86 

.76 

.76 

*81 

*85 

.425 

.85 

.75 

.78 

.87 

.73 

*72 

*81 

*79 

.475 

.82 

.71 

.83 

.93 

.82 

♦  75 

*89 

*76 

.525 

,93 

.73 

.86 

.94 

.74 

•  80 

*91 

*72 

.575 

.91 

.  82 

.91 

.84 

.71 

.73 

♦  95 

*80 

.625 

.90 

.84 

.87 

.79 

.83 

.75 

♦  87 

♦  71 

.675 

.93 

.  87 

.91 

.79 

.74 

.79 

*87 

*82 

.725 

.97 

.86 

.83 

.77 

.74 

.79 

*84 

*79 

.775 

1  .06 

.88 

.85 

.80 

.77 

.96 

.79 

♦  76 

.825 

1  .09 

.90 

.80 

.77 

.74 

.90 

.70 

♦  72 

.875 

.93 

.91 

.72 

.70 

.77 

.88 

♦  77 

.81 

.925 

.85 

.  85 

.79 

.75 

.78 

.91 

.82 

*74 

.975 

.76 

.79 

.84 

.77 

.82 

.90 

.81 

*81 
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TABLE  6 

ATTENUATION  factors  for  e0=  75  degrees,  the  fractional 
STANDARD  DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  18  PERCENT 
FOR  THE  ATTENUATION  FACTORS  DUE  TO  SCATTERED  RADIATION  ONLY* 
DATA  IN  THE  TABLES  SHOULD  BE  MULTIPLIED  BY  O.'  TO  GET  THE 
CORRECT  ATTENUATION  FACTOR. 


.49 

* 

.98 

1.47 

1.96 

2.94 

3.92 

4.90 

5*88 

I 

SCATTERED 

RADIATION 

.025 

.08 

.16 

.21 

.  13 

.20 

.17 

♦  18 

•  16 

.075 

.14 

.  16 

.15 

.  14 

.21 

.18 

.  19 

.  16 

.125 

.07 

.13 

.20 

.  15 

.  15 

.21 

.  19 

*  17 

.  175 

.12 

.  13 

.20 

.  17 

.20 

.13 

.  15 

*  16 

.225 

.  14 

.  17 

.13 

.21 

.20 

.19 

*  14 

.  18 

.275 

.15 

.  15 

.17 

.20 

.22 

.14 

♦  22 

.17 

.325 

.19 

.24 

.27 

.  14 

.  13 

.18 

*  IS 

.  16 

.375 

.24 

.24 

.17 

.19 

•  17 

.18 

.16 

.22 

.425 

.18 

.  18 

.  18 

.  17 

.  16 

.14 

.14 

•  20 

.475 

.16 

.  17 

.  18 

.  15 

.  14 

.12 

.20 

*  14 

.525 

.22 

.  17 

.14 

.  18 

.  15 

.19 

.16 

.23 

.575 

.28 

.16 

.  14 

.  19 

.20 

.16 

.20 

.20 

.625 

.15 

.24 

.17 

.17 

.  17 

.19 

.  19 

.18 

.675 

.30 

.15 

.16 

.  19 

.19 

.17 

.16 

.19 

.725 

•  21 

.13 

•  18 

.21 

•  19 

.16 

.16 

.19 

.775 

.27 

.  17 

.20 

.17 

.14 

.2! 

♦  14 

•  1  1 

.825 

.17 

.  19 

.21 

.  15 

.13 

•  20 

.18 

*21 

.875 

.19 

.20 

.  14 

.16 

.  16 

.18 

.  1  7 

*  17 

.925 

.18 

.  18 

.14 

.16 

.20 

*21 

.  18 

.15 

.975 

.08 

.20 

.17 

.28 

.22 

.21 

.21 

118 

SCATTERED 

PLUS  UNSCATTERED  RADIATION 

.025 

.09 

.16 

.22 

.  14 

.20 

.18 

.18 

416 

.075 

.14 

.16 

.15 

•  15 

•  21 

.18 

.  19 

.  1  6 

.  125 

.07 

.  14 

.21 

.  15 

«  15 

.21 

*  19 

i  18 

.175 

.12 

.  1  3 

•  20 

.  17 

.20 

.13 

*  15 

K  *7 

.225 

.14 

•  1  7 

.13 

.21 

*20 

.19 

i  15 

•  18 

.275 

.16 

.15 

.17 

.21 

.22 

.15 

.23 

1 17 

.325 

.19 

.24 

.28 

.  15 

♦  14 

.  18 

.20 

tie 

.375 

.24 

.24 

.17 

.19 

.  17 

*  19 

.  17 

ika 

.425 

•  18 

.  18 

.18 

.  18 

.  17 

.15 

*  14 

.20 

.475 

.16 

.  17 

.18 

.15 

.  14 

.12 

.20 

♦  15 

.525 

.22 

.  17 

.15 

.  19 

.  15 

.20 

.  17 

.23 

.575 

.28 

O  17 

.  14 

.  19 

.20 

.17 

•  20 

.20 

.625 

.16 

.24 

.17 

.  17 

.  17 

.19 

.  19 

.19 

.675 

.30 

.  16 

.17 

.19 

.20 

.17 

.16 

.19 

.725 

.22 

.13 

.19 

.22 

.20 

.17 

.16 

.  19 

.775 

.28 

.  18 

,21 

.  17 

.  14 

.21 

.14 

.  1  1 

.825 

.17 

•  20 

.21 

.15 

.  14 

.20 

.19 

.21 

.875 

.20 

.20 

.  15 

.  16 

.  17 

.18 

.18 

.  18 

.925 

.18 

.  19 

.  14 

.  17 

.20 

.21 

.  19 

.  16 

.975 

.09 

.21 

.17 

,28 

.22 

.21 

.21 

.  18 
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TABLE  7 

Comparison  of  ribbed  slab  and  equivalent  plane  slab  attenuation 
factors  for  scattered  radiation. 


mm 

Attenuation  Factors 

0 

Ribbed 

Plane 

Ratio 

Slab 

Slab 

0° 

0. 184±0. 001 

0. 198±0. 005 

0.93  ±  0.02 

45° 

0. 121±0. 003 

0. 109±0. 003 

1.11  ±  0.04 

o 

o 

tO 

0.  0640±0. 0007 

0. 0559±0.  0020 

1.14  ±  0.04 

75° 

0. 0176±0. 0007 

0.  0123±0. 0006 

1.43  ±  0.09 

TABLE  8 

Equivalent  plane  slab  attenuation  factors 
for  uncollided  radiation. 


9 

Attenuation 

o 

Factor 

0° 

0.124 

45° 

0.0528 

o 

o 

tO 

0.0158 

75° 

0.000350 
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TABLE  9 

ATTENUATION  FACTORS  FOR  RADIATION  EMITTED  INTO  DIRECTIONS 
WITHIN  2.5  DEGREES  OF  0O=  45  DEGREES.  The  FRACTIONAL  STANDARD 
DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  8  PERCENT  FOR  THE 
ATTENUATION  FACTORS  DUE  TO  SCATTERED  RADIATION  ONLY.  DATA  IN 
THE  TABLES  SHOULD  BE  MULTIPLIED  BY  0.1  TO  GET  THE  CORRECT 
ATTENUATION  FACTORS. 


98  1.47  1.96  2.94 


SCATTERED  RADIATION 


3.92 


4.90  5.88 


.025 
.075 
.125 
.  175 
.225 
.275 
.325 
.375 
.425 
.475 
.525 
.575 
.625 
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1 .90 


15 


done  for  the  equivalent  plane  slab  using  the  same  computer  code  that  was 
used  for  the  ribbed  slab  calculation.  10,000  histories  were  processed 
for  each  angle  of  incidence.  The  attenuation  factor  for  the  equivalent 
plane  slab  and  the  average  attenuation  factor  for  the  ribbed  slab  are 
compared  in  Table  7.  The  results  indicate  that  for  normal  incidence 
the  attenuation  factor  for  scattered  radiation  is  about  7%  less  for  the 
ribbed  slab  than  for  the  equivalent  plane  slab.  Or  the  other  hand,  for 
8q  =  75°,  the  ribbed-slab  attenuation  factor  is  about  40%  greater  than 
the  plane-slab  attenuation  factor.  For  the  sake  of  completeness,  the 
equivalent  plane  slab  attenuation  factors  for  the  uncollided  radiation 
are  presented  in  Table  8. 


Results  for  an  Imperfectly  Collimated  Plane  Source 


The  gamma  ray  source  used  in  the  experimental  ribbed-slab  studies 

did  not  produce  a  perfectly  collimated  beam  of  radiation.  In  order  to 

approximate  the  effect  of  the  diverging  beam  on  the  attenuation  factor, 

consider  a  plane  source  which  emits  radiation  isotropically  for  42.5°  £ 
o  6  6 

0  £  47.5  and  -  —  £  co  ^  -r  .  If  the  ribbed  slab  is  not  present,  then 

O  A  O  A 

the  exposure  measured  by  a  detector  near  the  plane  source  will  be  pro¬ 
portional  to 


o  cos9 


,  ,  ,  cos42.5 

dcp  =  6  In  - 

o  ro 

cos47. 5 


With  the  same  proportionality  constant,  the  exposure  when  the  ribbed 


slab  is  present  will  be  proportional  to 

f ',7-5°  f  1  a  <e  H,X) 

d9  sin9  dtp  — — — 


d9  sin8 
o  o 


Af(8Q,H,X) 


The  attenuation  factor  for  the  diverging  beam  can  then  be  defined  as 


A.(A9  ,H,X)  = -  d9  sin9 

f  °  °  ° 


Af(9o,H,X) 


cos42,5 

cos47.5< 


To  evaluate  the  integral,  additional  attenuation  factors  were  computed 

for  radiation  incident  on  the  ribbed  slab  with  9  =  42.5°  and  9  =  47.5°. 

o  o 
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For  each  of  these  angles  10,000  histories  were  analyzed.  Then,  using 
attenuation  factors  for  0q  =  42.5°,  45°,  and  47.5°,  the  integral  was 
evaluated  numerically  using  the  trapizoidal  rule.  Discontinuities  in 
the  derivative  of  the  integrand  were  ignored  in  the  calculation. 

Data  for  A^(A0q,H,X)  are  given  in  Table  9  and  are  presented  graphically, 
in  Figure  7.  Irregularities  due  to  uncollided  radiation  diminish  with  in¬ 
creasing  H,  in  contrast  to  the  situation  for  a  perfectly  collimated  beam. 

The  rate  with  which  ^.y(A0o,H,X)  loses  its  dependence  on  X  as  H  increases  is 
expected  to  be  greater,  the  larger  the  angle  of  divergence  of  the  beam  of 
incident  radiation. 

IV.  CONCLUSION 

When  the  height  of  the  detector  above  the  slab  is  about  equal  to  the 
rib  height,  the  data  in  this  report  show  that  the  attenuation  factor  for 
scattered  radiation  depends  strongly  on  the  horizontal  detector  position. 
However,  when  the  detector  height  above  the  end  of  the  ribs  becomes  equal 
to  or  greater  than  an  appreciable  fraction  of  the  rib  separation  distance, 
the  attenuation  factor  appears  to  be  independent  of  the  horizontal  posi¬ 
tion.  It  is  of  interest  to  note  that  if  the  attenuation  factor  is  aver¬ 
aged  over  one  rib-separation  distance,  the  average  value  is  not  equivalent 
to  the  corresponding  attenuation  factor  for  a  plane  slab  having  the  same 
average  mass  thickness  as  the  ribbed  slab.  For  radiation  normally  inci¬ 
dent  on  the  slab  the  average  attenuation  factor  for  the  ribbed  slab  is 
about  7%  smaller  than  that  for  the  equivalent  plane  slab.  On  the  other 
hand,  for  radiation  incident  at  grazing  angles,  the  data  indicate  that 
the  average  attenuation  factor  for  the  ribbed  slab  is  about  40%  greater 
than  the  value  for  the  equivalent  plane  slab. 

When  uncollided  radiation  is  included  in  the  attenuation  factor, 
then  it  is  found  that  for  a  perfectly  collimated  plane  parallel  beam 
of  radiation  incident  on  the  slab,  the  attenuation  factor  retains  a 
dependence  on  the  horizontal  detector  position  for  all  detector-slab 
separation  distances.  When  the  beam  of  radiation  incident  on  the  slab 
contains  radiation  traveling  in  a  small  cone  of  directions,  then  the 
total  attenuation  factor  slowly  loses  its  dependence  on  the  horizontal 
detector  position  as  the  detector-slab  separation  distance  increases. 
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Since  the  data  in  this  report  are  for  a  particular  ribbed  slab 
configuration,  it  is  not  possible  to  make  broad  generalizations  about 
the  effects  of  inhomogeneities  in  shields.  However,  the  present  re¬ 
sults  indicate  that  in  sane  cases  there  are  sufficiently  large  differ¬ 


ences  between  the  shielding  capabilities  of  inhomogeneous  slabs  and 
homogeneous  slabs  with  equivalent  average  mass  thicknesses  to  make 
further  study  desirable. 
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A  cross-sectional  view  of  the  ribbed-slab 
defining  the  various  parameters  used  to 
describe  the  geometry. 


NUMBER  PER  0.1  MeV 


FIGURE  2 

Energy  spectrum  of  radiation  assumed  incident  on  the  ribbed  slab, 


FIGURE  3 


Attenuation  factors  for  scattered  radiation  when  0 


Attenuation  factors 
included. 


FIGURE  d 


Attenuation  factors  for  0  =  60°  with  uncollided  radiation 

included. 


FIGURE  7 


Attenuation  factors  for  a  beam  of  radiation  incident  on  the 
ribbed  slab  with  directions  of  incidence  diverging  2.5°  on 
either  side  of  0O  =  45°.  Uncollided  radiation  is  included. 
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APPENDIX  A 

Calculation  for  a  Simulated  Wood  Floor 

The  wood  was  assumed  to  have  photon  interaction  properties  of  con- 

3 

crete  having  a  density  of  0.641  g/cm  .  The  linear  dimensions  were  chosen 
to  correspond  closely  to  a  1-1/2  inch  floor  supported  by  2-  by  10-inch 
joists  spaced  16  inches  apart  (measured  center  to  center).  These  dimen¬ 
sions  require  that 

x  =  8  inches; 

r 

x  =  7  inches; 

s 

z  =  1-1/2  inches; 

s 

z  -  11-1/2  inches, 

r 

The  source  radiation  was  assumed  to  have  the  energy  1.25  MeV  and  to  be 
incident  normally  on  the  floor. 

The  data  for  this  calculation  are  presented  in  Table  Al.  For  this 
case,  1 00,000  photon  histories  were  analyzed.  The  calculation  required 
less  than  four  minutes  on  the  IBM-7094.  Like  the  attenuation  factors 
for  scattered  radiation  presented  in  Section  III,  the  attenuation  factors 
for  the  simulated  wood  floor  exhibit  a  strong  dependence  on  horizontal 
detector  position  when  the  detector  height  is  about  equal  to  the  rib 
height.  When  the  detector  height  above  the  ribs  exceeds  an  appreciable 
fraction  of  the  rib  separation  distance,  the  dependence  on  the  hori¬ 
zontal  detector  position  is  no  longer  evident.  Unlike  the  attenuation 
factors  for  scattered  radiation  given  in  Section  III,  for  detector  heights 
about  equal  to  the  rib  height,  the  attenuation  factors  for  the  simulated 
wood  floor  have  their  maximum  value  for  detector  positions  over  the  ribs. 
This  results  because  the  floor  thickness  is  on  the  order  of  0.1  mfp  while 
the  rib  heights  is  on  the  order  of  1  mfp.  Thus,  radiation  incident  under 
the  ribs  is  more  likely  to  be  scattered  before  leaving  the  ribbed  slab. 
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TABLE  A 1 

ATTENUATION  FACTORS  FOR  1.25  MEV  GAMMA  RADIATION  INCIDENT 
NORMALLY  ON  A  SIMULATED  WOOD  FLOOR.  THE  FRACTIONAL  STANDARD  1 
DEVIATION  HAS  AN  AVERAGE  VALUE  OF  ABOUT  5  P5RCENT  FOR1  THE 
ATTENUATION  FACTORS  DUE  TO  SCATTERED  RADIATION  ONLY.  DATA  IN 
THE  TABLES  SHOULD  BE  MULTIPLIED  BY  0.1  TO  GET  THE  CORRECt 
ATTENUATION  FACTOR.  ' 
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.62 

.78 

.91 

1.16 

1  .41 

1.91 

2.16 

2*91 

SCATTERED  RADIATION 

.025 
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1.17 
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1.31 
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!  *  15 
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1.18 
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1.28 

1.17 
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1.27 

1 .19 
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1*26 
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1.21 
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1.36 
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1 .24 

1  .24 

1  .29 
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.325 
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1.23 
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1  .26 

1.30 

1.34 

.375 

1.00 

1.18 
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1.25 
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1  .31 

1.29 
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1.19 
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.475 

1.18 

1.19 
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1.28 

1.28 
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1.17 

1.19 
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1.19 

1.14 
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1  <,28 

1  .33 

1  «  19 
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1.21 

1  .23 

1.30 
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1  .28 

1.25 
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1.15 
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1.31 
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1.16 

1.15 
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1.19 
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1.20 

1.29 
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1.13 

1.32 

1.24 
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1  .  18 

1  .26 
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1  *  19 
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1 .19 
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1  .27 
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1 .31 

1  .26 
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1  .31 

1  .26 

1  .  3® 

1.25 

1.27 

1  .28 

1.30 

1*16 

.875 

1  .43 

1  .37 

1.27 

1.27 

1 .29 

1  .37 

1*32 

1*21 

.925 

1  .51 

1.38 

1.34 

1.25 

1 .28 

1  .29 

1*25 

1  *27 

.975 

1  .72 

1  .39 

1.29 
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1  .30 

1  *21 

1  .23 

1 124 

j  SCATTERED 

PLUS  unscattered  RADIATION 

'  1  1 

.025 

5.00 

4.89 

4.83 

4.62 

4.70 

4*76 

4*75 

4  1  68 

.075 
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9.85 

9.87 

1  ©*03 
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.  125 

iO.Ol 

9.93 

9.88 
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9.97 

9.91 

1  0.00 
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9.86 

9.99 

9.99 

9,87 

9.99 

9,95 

9*93 

9198 

.225 

9.87 

9.92 

9.90 
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9.90 
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9.96 
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9.91 

10.07 

9.89 

9.94 
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9.99 

9191 

.325 

9.90 

9.77 

9.93 

9.94 

9.96 

9.96 
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.375 

9.71 

9.89 

9.92 

9.95 
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9.87 

9.98 

9.92 

.425 

9.87 

9.92 

9.90 

10.01 

9.99 

9.95 

9.90 

9.97 

.475 

9.89 

9.90 

9.91 

9.98 

9.98 
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9.87 

9.89 

.525 

9.89 
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9.92 
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9.98 
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9.89 

10.04 

.575 

9.79 

9.91 

9.93 
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9.95 

9.90 

9.98 

9.96 

.625 

9.86 

9.85 

9.94 

9,94 

9.90 

9.92 
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10.02 

.675 

9.86 

9.85 

9.88 

9.89 

9.98 

9.97 

9.91 

10*09 

.725 

9.83 

10.03 

9.94 

9.99 

9.88 

9.96 

9.89 

9.90 

.775 
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10.03 
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9.92 
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9.97 

9.95 

9.99 

.825 
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9.95 

9.97 

9.98 

10.00 

9.87 

.875 

10.14 
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9.98 
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9.99 

10.07 
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9.92 

.925 

10.22 
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9.98 

10.00 

9.95 

9.97 
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5.17 

4.84 

4.74 

4.76 

4.75 

4 . 66 

4*68 

4*69 
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APPENDIX  B 


1.  Description  of  the  Computer  Program. 

A  general  description  of  the  calculation  of  the  data  in  this  report 
has  been  given  in  Section  II.  A  somewhat  more  detailed  description  of 
the  individual  subroutines  which  make  up  the  computer  program  will  now 
be  undertaken. 

ADAM  RIB  II  (Main  Program) 

Subroutines  called:  LEARN 3 ,  DATEX3 ,  TRACK3,  TEACH3. 

The  main  program  consists  of  three  loops.  The  inner  loop  which 
the  program  executes  IHPP  times  for  each  value  of  the  pair  L  and  NXS0 
determines  the  number  of  histories  to  be  done  for  each  source  point. 

NXS0  determines  the  location  of  a  source  point  and  NXSX  is  the  maxi¬ 
mum  number  of  source  points.  L  is  an  index  which  assumes  all  values 
from  1  to  L00P.  L00P  is  defined  below.  The  product  L00P*  NXSX*  IHPP 
determines  the  total  number  of  histories  to  be  processed. 

Subroutine  LEARN3 

Subroutines  called:  none. 

This  subroutine  reads  all  input  data  and  prints  them  under  appro¬ 
priate  labels.  The  first  card  read  contains  72  alphanumeric  characters 
which  when  printed  identify  the  run.  The  input  variables  will  be  de¬ 
fined  in  the  order  in  which  they  appear  in  the  subroutine. 

IHIX:  The  number  of  histories  to  be  processed  for  a  given  value  of 

the  index  L  in  the  main  program. 

NPAX:  The  number  of  detector  heights  to  be  considered  in  the  cal¬ 

culation. 

NDAX:  The  number  of  detectors  to  be  considered  at  a  given  height. 

NMAX:  The  maximum  number  of  interactions  allowed  in  a  given  history. 

For  the  calculations  in  this  report,  this  parameter  was  assign¬ 
ed  the  value  50, 

NMIN:  The  minimum  number  of  interactions  which  a  photon  had  to  have 

before  it  was  allowed  to  contribute  to  the  detector  response. 
For  the  calculations  in  this  report,  this  parameter  was  assign¬ 
ed  the  value  1.  This  meant  that  the  Monte  Carlo  calculation 
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MX  AX: 

NXSX: 


NESX: 

L00P: 


NRAN: 


XS: 

XRJ 

ZS: 

ZR: 

ZP: 

UZCT: 

CTH0: 

PH0: 


DEN: 


EB: 


DUM1 

DUM2 

DUM3 


gave  the  detector  response  due  to  scattered  radiation  only. 
The  number  of  energies  for  which  interaction  coefficients 
are  read. 

The  number  of  source  points  to  be  located  between  the  cen¬ 
ters  of  two  ribs.  The  distance  between  rib  centers  was 
divided  into  NXSX  intervals  and  a  source  point  was  sel¬ 
ected  from  the  center  of  each  interval. 

The  number  of  source  energies  which  makes  up  the  source 
spectrum. 

The  number  of  times  the  main  program  executes  its  outer 
loop.  L00P*  IHIX  is  the  total  number  of  histories  to  be 
processed  in  the  run. 

T'j  initial  random  number  read  in  an  octal  format,  Rules 
for  selecting  this  number  are  given  in  the  program  list- 
'ing  for  the  random  number  subroutine  (RAND,  A,B,C,D). 

Half  the  distance  between  adjacent  edges  of  two  adjacent 
ribs. 

Half  the  distance  between  the  centers  of  two  adjacent  ribs. 
The  thickness  of  the  slab. 

The  thickness  of  the  slab  plus  the  height  of  a  rib. 

A  list  of  detector  heights  measured  from  the  side  of  the 
slab  opposite  the  ribs. 

The  cosine  cutoff  described  in  Section  II,  B. 

The  initial  value  of  the  direction  cosine  u  . 

z 

This  angle,  read  in  degrees,  defines  the  initial  direction 
cosine  relative  to  the  x-axis.  We  have 

u  =  ~\j  1  -  CHT02  cos  (PH0) . 

The  density  of  the  ribs  and  the  slab.  In  the  calculation 
described  in  this  report,  the  units  for  this  quantity  were 
grams  per  square  centimeter  per  inch. 

A  list  of  energies  in  MeV  for  which  interaction  coefficients 
are  to  be  read. 


These  variables  are  used  to  skip  interaction  coefficient  data 
on  the  input  cards  which  are  not  used  in  the  calculation. 
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XSECB:  The  array  containing  the  input  interaction  coefficient  data. 

The  pair  of  indices  (M,l)  refers  to  the  Compton  interaction 
coefficient,  (M,2)  to  the  pair  production  interaction  coeffi¬ 
cient,  (M,3)  to  the  total  attenuation  coefficient,  and  (M,4) 
to  the  energy  absorption  coefficient  for  air.  For  the  calcu¬ 
lations  described  in  this  report,  all  these  quantities  had 
units  square  centimeter  per  gram. 

EMIN:  The  cutoff  energy. 

ES0:  A  list  of  energies  representing  the  source  spectrum. 

WTS0:  A  list  of  probabilities  assigned  to  the  list  of  source  ener¬ 

gies  ES0. 

Subroutine  DATEX3 

Subroutines  called:  TABIN 

This  subroutine  performs  a  number  of  preliminary  tasks  in  prepar¬ 
ation  for  the  main  part  of  the  calculation.  First,  the  interaction 
coefficient  data  are  modified  so  that  they  will  contain  data  more 
directly  useful  in  later  stages  of  the  calculation.  XSECB(M,1)  is 
changed  to  the  probability  that  given  an  interaction,  either  a  Comp¬ 
ton  scattering  or  a  pair  production  interaction  will  occur.  XSECB 
(M,2)  is  multiplied  by  the  density  DEN  to  convert  it  to  units  con¬ 
sistent  with  the  units  used  in  specifying  the  dimensions  of  the  slab 
and  ribs. 


Second,  the  logarithm  of  each  entry  in  the  array  XSECB  and  the 
list  EB  is  calculated.  Double  logarithmic  interpolation  is  then 
used  to  expand  the  tables  of  interaction  probabilities.  Entries  in 
these  expanded  tables  correspond  to  energies  defined  by  the  formula 

E  _  — 1000 -  _  i  .  1  <;  NE  £  899 

NE  +  90.5 


In  subsequent  parts  of  the  program,  if  an  interaction  probability 
is  desired  for  a  certain  photon  energy  E  ,  the  index  NE  is  defined 
by  rounding  the  following  point  number 


1000 
E  +  1 


90 


down  to  the  nearest  whole  number.  Then  the  interaction  probability 
corresponding  to  the  index  NE  is  taken  as  the  interaction  probability 
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for  photons  of  energy  E  .  This  particular  table-look-up  algorithm 

4 

was  suggested  by  A.  B.  Chilton  . 

Third,  a  table  of  sines  and  cosines  is  generated  for  360  angles 
o  o 

starting  at  0.5  and  proceeding  in  steps  of  1  up  to  and  including 

o  o  o 

359.5  .  The  sine  and  cosine  of  a  randon  angle  between  0  and  360 

are  then  selected  by  computing  an  index  according  to  the  formula 
NPH  =  360.0*RAN,  where  RAN  is  a  random  number  uniformly  distribu¬ 
ted  on  the  interval  (0,1)  and  then  choosing  the  sine  and  cosine 
in  this  table  corresponding  to  the  index  NPH  +  1  . 

Fourth,  the  quantities  XR2  and  XS2  are  computed.  These 
give  the  distance  between  rib  centers  and  the  separation  distance 
between  adjacent  edges  of  the  ribs,  respectively. 

Fifth,  the  cutoff  energy  EMIN  is  converted  to  its  corres¬ 
ponding  Compton  wavelength  WMAX  . 

Sixth,  a  list  of  source-point  coordinates  XS0  is  computed. 
These  are  defined  by  dividing  the  interval  between  the  centers  of 
two  adjacent  ribs  into  NXSX  increments  and  selecting  the  midpoint 
of  each  increment  as  a  source  point. 

Seventh,  a  check  is  made  to  see  if  the  input  parameter  IHIX 
is  divisible  by  the  number  of  source  points  NXSX  .  If  this  is  not 
the  case,  then  the  program  terminates. 

Eighth,  IRC  is  set  equal  to  the  initial  random  number  IRAN  , 
the  sine  of  the  angle  of  incidence  is  computed,  and  the  sine  and  co¬ 
sine  of  the  angle  PH0  are  calculated. 

Ninth,  the  arrays  and  lists  which  are  used  “o  accumulate  in¬ 
dividual  scores  and  the  squares  of  individual  scores  are  initially 
set  equal  to  zero.  The  names  of  these  variables  are  defined  as 
follows: 

TS2N:  This  variable  is  used  to  accumulate  the  statistical  weight 

of  a  photon  when  it  is  transmitted.  After  division  by  the 
total  r umber  of  histories,  it  gives  the  number  current  of 
photons  transmitted  by  the  ribbed  slab. 

TF2N:  This  variable  is  used  to  accumulate  individual  contributions 

to  the  transmitted  number  flux. 

TFEX:  This  variable  is  used  to  accumulate  individual  contributions 

to  the  exposure  due  to  transmitted  photons. 
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BSCN:  This  list  is  used  to  accumulate  contributions  to  the  reflect¬ 

ed  number  current  as  a  function  of  the  position  of  emergence 
relative  to  the  ribs. 

BSFN:  This  list  is  used  to  accumulate  contributions  to  the  reflect¬ 

ed  number  flux  as  a  function  of  the  position  of  emergence. 
EBSC:  This  list  is  used  to  accumulate  contributions  to  the  exposure 

made  by  reflected  photons,  again  as  a  function  of  the  point 
of  emergence. 

TS1N:  This  array  is  used  to  accumulate  contributions  to  the  trans¬ 

mitted  number  current  as  a  function  of  the  height  of  the 
detector  above  the  slab  and  the  horizontal  position  of  the 
detector  relative  to  the  ribs. 

TSFN:  This  array  is  used  to  accumulate  contributions  to  the  ex¬ 

posure  due  to  transmitted  photons  as  a  function  of  detector 
height  and  horizontal  detector  position. 

The  cosine  cutoff  described  in  Section  II  is  applied  only  to 
the  arrays  TSFN  and  TESC  .  For  each  of  the  variables  listed  above 
there  is  a  second  variable  named  by  attaching  the  symbol  2  to  the 
end  of  the  names  listed.  These  variables  are  used  to  accumulate  the 
squares  of  individual  contributions  so  that  standard  deviations  can 
be  calculated  at  the  end  of  the  run. 

Tenth,  the  probabilities  for  individual  source  energies  are 
used  to  calculate  a  cumulative  probability  distribution  for  the 
source  energies.  This  cumulative  distribution  is  then  used  in  the 
random  sampling  of  the  initial  photon  energy.  The  program  normalizes 
the  probability  distribution  and  the  resulting  cumulative  distribution 

Subroutine  TRACK3 

Subroutines  called:  RANDC,  RANDA,  CHECK3 ,  C0MPT3 ,  GRADE3. 

This  subroutine  generates  path  lengths  between  interaction  points 
and  decides  what  kind  of  interactions  take  place.  The  first  part  of 
the  program  establishes  the  initial  photon  energy,  direction,  and  posi 
tion.  The  random  numbers  used  within  a  given  history  are  generated 
by  subroutine  RANDA  .  However,  this  chain  of  random  numbers  is  init¬ 
ialized  at  the  beginning  of  each  history  using  a  random  number  gener¬ 
ated  by  subroutine  RANDC.  The  total  number  of  random  numbers  gener¬ 
ated  by  RANDC  will  be  equal  to  the  number  of  histories.  The  y- 
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coordinate  is  defined  and  carried  along  in  this  subroutine,  but  this 
is  unnecessary  and  could  be  eliminated.  WT  is  the  statistical  weight 
of  the  photon  and  is  initially  set  equal  to  unity.  The  integer  N  is 
used  to  count  the  number  of  interactions  in  a  given  history.  When  a 
photon  escapes  from  the  slab,  N  is  compared  with  NMIN  to  see  if 
a  score  should  be  recorded.  If  N  exceeds  NMAX  ,  the  history  is 
terminated. 

Once  the  initial  conditions  have  been  defined,  the  track 
length  to  the  next  interaction  is  sampled  and  the  coordinates  of 
the  location  of  this  interaction  are  computed.  Then  subroutine 
CHECK3  is  called  to  determine  if  any  boundaries  have  been  crossed. 

If  a  boundary  has  been  crossed  CHECK3  makes  appropriate  adjust¬ 
ments  on  the  location  of  the  interaction  and  defines  NSCT  .  If 
NSCT  is  negative;  then  the  photon  has  escaped  and  TRACK3  decides 
whether  or  not  a  score  should  be  recorded.  If  NSCT  =  0,  then 
CHECK3  found  that  the  position  of  the  interaction  point  was  not 
in  the  ribbed  slab  but  that  when  the  photon  was  moved  along  the 
direction  it  was  going,  it  subsequently  re-entered  the  slab.  In 
this  case,  TRACK3  samples  a  new  track  length  to  move  the  photon 
away  from  the  boundary  of  the  ribbed  slab  and  into  the  slab  or 
rib. 

If  NSCT  >  0,  then  an  interaction  point  within  the  slab  has  re¬ 
sulted.  If  the  index  NE  is  equal  to  404,  the  photon  energy  is 
very  nearly  equal  to  the  threshold  energy  for  pair  production.  Thus, 
if  NE  <  404,  TRACK3  must  choose  between  a  pair-production  inter¬ 
action  and  a  Compton  scattering.  If  NE  >  404,  only  Compton  scatter¬ 
ing  is  allowed.  When  a  pair-production  interaction  is  selected, 
the  photon  energy  is  set  equal  to  0.511  MeV  (Compton  wavelength  is 
unity  and  NE  =  571),  its  statistical  weight  is  multiplied  by  two, 
and  a  new  direction  is  sampled  from  an  isotropic  angular  distribu¬ 
tion.  Sampling  a  new  direction  and  energy  in  the  case  of  Compton 
scattering  is  accomplished  by  subroutine  C0MPT3. 

Photoelectric  absorption  is  taken  into  account  by  multiplying 
the  statistical  weight  of  the  photon  by  SURV,  the  probability  that 
an  interaction  is  a  pair  production  or  a  Compton  scattering  given 
that  an  interaction  occurs. 
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Subroutine  CHECK3 

Subroutines  called:  none. 

This  subroutine  checks  to  see  if  boundaries  have  been  crossed. 
Most  of  the  questions  asked  by  the  subroutine  are  phrased  for  pho¬ 
tons  whose  direction  cosine  along  the  x-direction  is  positive.  If 
this  direction  cosine  is  negative,  then  the  transformation 

UX  -  -UX  =  PUX 

XN  -  -XN 

is  made  and  the  photon  is  treated  as  if  the  direction  cosine  were 
positive.  This  transformation  is  permissible  because  the  origin 
of  the  x-coordinate  is  located  halfway  between  the  ribs. 

Basically,  the  subroutine  checks  for  four  events.  First  it 
checks  to  see  if  the  photon  has  been  transmitted  by  the  ribbed 
slab.  If  it  has,  the  x-coordinate  of  the  point  where  the  photon 
crossed  the  plane  defined  by  the  top  of  the  ribs  is  computed  and 
NSCT  is  set  =  -1. 

Second,  the  subroutine  checks  to  see  if  the  photon  has  been 
reflected  by  the  slab.  If  reflection  has  occured,  the  x-coordin- 
ate  of  the  reflection  point  is  computed  dnd  NSCT  is  set  =  -1. 

Third,  the  subroutine  checks  to  see  if  the  photon  has  left 
the  ribbed  slab.  If  this  is  the  case,  then  either  the  photon 
has  actually  been  transmitted  or  the  photon  will  re-enter  the  rib¬ 
bed  slab.  If  transmission  has  occured,  then  the  x-coordinate  of 
the  point  where  the  photon  crosses  the  plane  defined  by  the  top  of 
the  ribs  is  computed  and  NSCT  is  set  =  -1.  Otherwise  the  coordin¬ 
ates  of  the  point  where  the  photon  re-enters  the  slab  are  determined 
and  NSCT  is  set  =  0. 

Fourth,  the  subroutine  checks  to  see  if  the  x-coordinate  of 
the  photon  has  exceeded  the  x-coordinate  of  the  rib  center.  If  this 
has  happened,  then  the  separation  distance  between  the  rib  centers 
is  subtracted  from  the  x-coordinate  of  the  photon  and  the  question 
is  asked  again.  This  process  is  continued  until  an  x-coordinate  is 
obtained  which  is  less  than  the  x-coordinate  of  the  rib  center. 
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Once  the  subroutine  determines  that  the  photon  has  not  been  re¬ 
flected  or  transmitted  by  the  slab,  and  that  the  photon  has  not  left 
and  re-entered  the  slab,  then  it  sets  NSCT  =  1  , 

Subroutine  C0MPT3 

Subroutines  called:  RANDA. 

5 

This  subroutine  uses  the  method  of  Kahn  to  sample  a  new  photon 
direction  and  energy  from  the  Klein-Nlshina  distribution. 

Subroutine  GRADE3 

Subroutines  called:  none. 

This  subroutine  computes  scores  in  the  event  that  a  photon  is  re¬ 
flected  or  transmitted  by  the  ribbed  slab.  The  variables  in  which  the 
scores  are  accumulated  are  defined  in  the  description  for  subroutine 
DATEX3. 

If  a  photon  is  reflected,  its  contribution  to  the  reflected  num¬ 
ber  current,  number  flux,  and  exposure  is  computed  as  a  function  of 
the  position  of  emergence  of  the  photon.  Because  the  albedo  was  of 
only  minor  interest  in  the  present  study,  the  additional  complication 
of  the  cosine  cutoff  discussed  in  Section  II  was  not  introduced  in  the 
computation  of  reflected  number  flux  and  exposure. 

When  a  photon  is  transmitted,  its  contribution  to  the  transmitted 
number  current,  number  flux,  and  exposure  is  computed  in  two  ways. 
First,  the  contribution  is  recorded  without  taking  into  account  the 
x-coordinate  of  the  photon  when  it  left  the  ribbed  slab.  In  this 
case,  the  cosine  cutoff  described  in  Section  II  was  not  used  in  the 
calculation  of  number  flux  and  exposure.  Second,  the  contribution 
is  recorded  as  a  function  of  detector  height  above  the  slab  and 
horizontal  detector  position.  In  this  second  case,  the  cosine  cutoff 
described  in  Section  II  was  used  in  the  calculation  of  number  flux 
and  exposure.  By  averaging  the  number  flux  and  exposure  (calculated 
using  the  second  method)  over  all  horizontal  detector  positions  and 
comparing  with  the  results  obtained  using  the  first  method  it  was 
possible  to  estimate  the  error  introduced  by  the  cosine  cutoff. 
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Subroutine  TEACH3 

Subroutines  called:  RIBTAB,  DIRIB. 

This  subroutine  normalizes  and  prints  the  output  data.  It  also 
computes  fractional  standard  deviations  for  all  results.  Headings  are 
printed  with  all  the  data  so  that  the  output  is  fairly  self-explana¬ 
tory. 

In  addition,  this  subroutine  calls  subroutine  DIRIB  which  com¬ 
putes  the  contribution  to  the  exposure  due  to  uncoil iriod  radiation. 
These  results  are  then  added  to  the  data  obtained  in  fche  Monte  Carlo 
calculation.  Because  of  an  oversight  there  i  i  no  provision  in  the 
program  for  determining  whether  the  direct  contribution  should  be 
added  in.  If  for  example  the  direct  radiation  is  included  in  the 
Monte  Carlo  calculation,  then  TEACH3  will  add  this  contribution  as 
calculated  by  DIRIB  anyway  with  the  result  that  some  of  the  output 
tables  will  be  incorrect. 

The  subroutine  also  punches  the  ribbed  slab  attenuation  factors 
and  the  errors  of  the  attenuation  factors  on  cards. 

Subroutine  DIRIB 

Subroutines  called:  TABIN 

This  subroutine  computes  the  contribution  to  the  exposure  due 
to  uncollided  radiation  as  a  function  of  the  horizontal  detector 
position  XD  and  the  distance  between  the  detector  and  the  source 
plane  ZP  .  The  result  is  returned  to  TEACH3  by  means  of  the 
argument  DIR  .  The  direction  of  the  uncollided  radiation  is  speci¬ 
fied  by  CPH0  and  CTH0  .  The  horizontal  detector  position  is 
taken  as  the  midpoint  of  the  detection  intervals  described  in  Section 
II. 

The  bulk  of  the  program  is  taken  up  with  the  calculation  of  the 
path  length  which  the  uncollided  radiation  must  travel  within  the 
ribbed  slab  in  order  to  reach  the  detector.  This  calculation  must 
take  account  of  the  fact  that  in  addition  to  passing  through  the 
basic  slab,  the  radiation  may  pass  through  parts  of  one  or  more  ribs. 
The  calculation  is  straightforward  but  rather  tedious. 
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After  the  path  length  T  has  been  computed,  the  last  part  of 
the  program  computes  the  contribution  to  the  exposure  made  by  each 
energy  in  the  source  spectrum. 

Subroutine  RIBTAB 

Subroutines  called:  none. 

This  subroutine  prints  tables  of  data  as  a  function  of  detector 
height  and  horizontal  detector  position.  Because  several  of  these 
tables  are  printed,  this  part  of  TEACH3  was  separated  into  a  sec¬ 
ond  subroutine  to  save  coding. 

The  maximum  table  width  is  eleven  columns.  If  more  than 
eleven  columns  are  needed,  then  the  additional  columns  are  includ¬ 
ed  in  a  second  table  which  is  printed  below  the  first  table. 

Subroutine  TABIN 

Subroutines  called:  none. 

This  subroutine  performs  a  quadratic  interpolation  simulta¬ 
neously  on  NMAX  functions  f  (x)  .  It  is  assumed  that  the  func- 

m 

tions  are  known  at  NMAX  values  x  .  The  program  assumes  that 

n 

the  values  x  are  stored  in  descending  order.  The  list  x  is 
n  n 

denoted  by  XB(N)  and  the  array  f  (x  )  by  FB(N,M)  ,  The  point 

m  n 

where  the  interpolated  values  are  to  be  found  is  given  by  X  and 
the  interpolated  values  by  F(M)  .  Unless  X  is  closest  to  the 
initial  or  final  entries  in  the  list  XB(N)  ,  the  central  point 
in  the  three-point  interpolation  is  taken  as  the  value  in  the  list 
which  is  closest  to  X  . 

The  first  time  the  program  is  called,  the  parameter  NTABIN 
should  have  the  value  one  .  Then  certain  differences,  products, 
and  sums  which  are  needed  each  time  the  subroutine  is  called  are 
computed  and  stored.  NTABIN  is  then  set  equal  to  two  and  on  sub¬ 
sequent  calls,  these  quantities  are  not  recomputed. 

Subroutine  RANDA,B,C,D 


Subroutines  called:  none 
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This  subroutine  is  actually  four  subroutines  in  one.  It  is 

used  to  calculate  pseudorandom  numbers,  uniformly  distributed  on 

the  interval  (0,1)  .  The  multiplicative  congruential  method 

35 

with  modulus  2  is  used.  If  the  calling  statement  is 
CALL  RANDA  (IR,R) 

15 

the  multiplier  is  5  .If  the  terminal  letter  on  the  sub- 

13  11  9 

routine  name  is  B,C,  or  D,  the  multiplier  is  5  ,5  ,  or  5  , 

respectively.  The  argument  IR  is  the  integral  form  of  the  ran¬ 
dom  number  and  R  is  the  normalized  form.  IR  must  be  supplied 
by  the  user  the  first  time  the  subroutine  is  called. 

Except  for  some  unimportant  comment  cards  at  the  beginning, 

this  subroutine  is  identical  to  one  described  in  more  detail  by 
0 

Spencer  . 

2.  Program  Listings  and  Sample  Output 

Listings  of  the  programs  described  above  are  given  on  the 
following  pages.  A  sample  output  from  the  program  follows  the 
listings. 
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FASTRAN 

ADAM  RIB  II  6-7-67 

COMMON  BSCN , BSCN2 ♦ BSFN , BSFN2  *  CPH .  CPHO . CTHO , DEN , EB , EBSC .  EBSC2  . 

1  EMIN.ESO.EXC.FNDAX. IHIX. IHPP, IRA, I RC , LOOP ,MXAX , NDAX , NE ♦ NESO, 

2  NESX . NMAX ,  NM I N , NPAX ♦ NRAN , NSCT , NXSO , NXSX , PCS , PHO ♦ RAN , S , SPH  » 

3  SPHO.STHO.SURV . SUZ ♦ TESC , TESC2 , TFEX » TFEX2 . TF2N , TF2N2 . TSFN , TSFN2 , 

4  TS 1 N . TS 1 N2 . TS2N , TS2N2 , UT , UX ♦ UY ♦ UZ , UZCT . W . WMAX . WT , WTSO , WTS0C . 

5  X  «  XN . XR , XR  S . XR2 . XS . XS2 , XSECB , XSO , Y , YN , Z ,  ZN , ZP  .  ZR . ZS 
DIMENSION  BSCN( 30) ,BSCN2(30 ) * BSFN (30 ) ,BSFN2(30 ) , CPH (360) . 

1  EB(25) .EBSC (30) ,EBSC2(30> ,ES0(30> ,EXC(899) .PCS (404) .SPH(360) ♦ 

2  SURV ( 899 ) «  TESC ( 30 . 30 ) . TESC2 ( 30 . 30 ) .TSFN (30 .30) , TSFN2 ( 30 ♦ 30 ) , 

3  TS1N(30.30)  .TSl N2 ( 30  *  30 ) »  UT  <  899  >  *WTS0(30). WTSOG ( 30 ) ♦ 

4  XSECB (25.4). XSO ( 80 ) , ZP ( 30 ) 

CALL  LEARN3 

CALL  DATEK3 
DO  30  L=l. LOOP 
NXSOaO 

10  NXSOaNXSO+1 

DO  20  J* 1 .IHPP 
20  CALL  TRACK3 

I F ( NXSO— NXSX ) 1 0 . 30 . 30 
30  CONTINUE 

CALL  TEAGH3 
CALL  SYSTEM 
END 
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fastran 

SUBROUTINE  LEARN3  6-12-67 

SUBROUTINE  LEARNS 

COMMON  BSCN, BSCN2 , BSFN . BSFN2  . CPH «  CPHO ♦  CTHO  4 DEN  4  EB  4  EBSC » EBSC2  4 
'  1  EMIN,ESO,EXC,FNDAX.  IHIX«  IHPP.  IRA»  IRC»L00P«MXAX»NDAX«NE,NES04 

2  NESX ♦ NMA X » NM I N , NPAX , NRAN«  NSC  I , NXSO . NXSX « PCS . PHO ♦ RAN . a , SPH . 

3  SPHO « STHO »  SURV  « SUZ  « TESC »  TESC2  »  TFEX  »TFEX2  *  TF2N «  TF2N2  « ISFN»TSFN24 

4  TSlN«TSl  N2 »  T  S2N «  TS2N2  « UT  *  UX  »UY  *  UZ  »  UZCT  *  W  «  wMAx  «  w  t  *  wToO»  w  i  oOC  « 

5  X ,  XN.  XR ,  XRS  ♦  XR2  .  XS  i  XS2  *  XSECB  *  XSO  »  Y  »  YN » Z  *  ZN «  ZP  •  ZR  *  ZS 
DIMENSION  BSCN(30) «BSCN2(30) iBSFNOO) »BSFN2(3u) «CPHi36u) , 

1  EB ( 25 ) ♦ EBSC (30) « EBSC2 ( 30 ) *  ESO( 30 )  »EXC(899) *  PCS (404) *SPH(360) » 

.  2  SURV (899) »  TESC ( 30  *  30 ) * TESC2 ( 30 ♦ 30 ) .  ISFN(3u»3u) ♦ 1 aFN2 i 30 . 3u ) « 

■  3  TS 1 N (30  4  30 ) *  TS 1 N2 ( 30  *  30 ) »  UT ( 899 ) »  WTSO( 30 ) «  WTSOC (30 ) 4 

4  XSECB ( 25  <4 )»XS0(80)»  Zp ( 30 ) 

WOT  6  » 10 
10  FORMAT ( 1 HO ) 

15  FORMAT ( 1  HI ) 

R I T7«  20 
WOT  6  «  20 
20  FORMAT (72H 

1  ) 

W0T6i 10 
W0T6  *  30 

30  FORMAT (54H  IHIX  NPAX  NDAX  NMAx  NM I N  MaAa  Naja  NE^a  LOOP) 
R IT7  »  40 1 I H I X  «  Np AX « NDAX  « NMAX  *  NM  JN • MX AX • NXSX  «  NESX • LOOP 
WOT  6  *  40 « IHIX* NPAX  » NDAX  «  NMAX *  NM 1 N  %  MxAx  »  Nxsx  « NESx • LOOP 
40  FORMAT! 1216) 

W0T6« 10 
W0T6 ♦ 50 

50  FORMAT (22H  INITIAL  RANDOM  NUMBER) 

RIT7i60,NRAN 
W0T6«60<  NRAN 
60  FORMAT (01 3) 

WOT 6  4  10 
W0T6  «  70 

70  FORMAT ( 31 H  XS  XR  ZS  ZR ) 

R I T7 « 80 1XS1 XR«  ZS  «  ZR 
WOT 6 »80»XS»XR4  ZS»ZR 
80  FORMAT (9F8«3 ) 

W0T6 « 10 
WOT  6  »  90 

90  FORMAT (61 H  DETECTOR  PLANE  HEIGHTS  (MEASURED  FROM  EN i RANCE  FACE  OF 
1  SLAB )  ) 

RIT7»80.  (ZP(  NP)  ,NPd  iNPAX) 

WOT618O, (ZP( NP)  ,NP=1 » NPAX ) 

WOT64 10 
WOT 6 » 100 

100  FORMAT  (  35H  COSINE  Cu  i  OFF  FOR  FLwx  COMP_  {"ATI  ON) 

RlT7»80«UZCT 
WOT 6  1  80  « UZC1 
WOT 6 » 10 
WOT 6  » 1 1 0 

110  FORMAT ( 29H  COSINE  OF  ANGLE  OF  INCIDENCE) 

R  I  T7 »  80 «  Cl  HO 
WOT 6  %  80 «  CTHO 
WOT61 1 0 
W0T6. 120 
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120  FORMAT { 25H  AZIMUTH  OF  INCIDENT  BEAM) 

R1T7»80«PH0 
W0T6  *80  «  PHO 
W0T6* 10 
W0T6* 130 

130  FORMAT (25H  DENSITY  OF  SLAB  AND  RIBS) 

RIT7«80*DEN 
W0T6  *80  *  DEN 
W0T6 ♦ l 5 
W0T6* 10 
RIT7* 140 
W0T6* 140 

140  FORMAT (80H 

1  ) 

W0T6* 10 
W0T6* 150 

150  FORMAT (37H  ENERGY  COMPTON  PAIR  TOTAL) 

DO  160  M= 1 « MXAX 

R I T7  « 1 70  *  EB  C  M ) , DUM 1 , XSECB  C  M , 1 ) , DUM2 ♦ XSECB  C  M ♦ 2 )  * DUM3 , XSECB ( M  *  3 ) 

1 60  WOT 6  *  1 70 . EB  <  M ) *  ( XSECB (M»l),l=l*3) 

170  FORMAT ( F7  *3*1  P7E 1 0  «  3 ) 

W0T6* 15 
WOT  6. 10 
WOT  6  « 1 80 

1B0  FORMAT (39H  ENERGY  ABSORPTION  COEFFICIENTS  FOR  AIR) 

WOT6* 10 

DO  190  M=1.MXAX 
RIT7,170»EB(M) «XSECB<M*4) 

190  WOT 6 *  170 • EB ( M) *  XSECB ( M  * 4 ) 

W0T6* 10 
WOT  6  «  200 

200  FORMAT ( 1 4H  CUTOFF  ENERGY) 

R I T7  «80  «  EM  I N 
W0T6»80«EMIN 
W0T6* 10 
WOT  6*210 

210  FORMAT (16H  SOURCE  ENERGIES) 

RIT7*80*  (ESO(NESO) *NES0=1 *NESX) 

W0T' *80* (ESO( NESO) ,NES0=1 *NESX) 

W0T6* 10 
W0T6*  220 

220  FORMAT (22H  ''URCE  ENERGY  ..EIGHTS) 

R I T7« 230  « ( WTSO( NESO) »NES0=1 .NESX) 

WOT6«230* (WTSO(NESO) *NES0=1 *NESX) 

230  FORMAT ( 1 P7E 1 0 • 2 ) 

WOT  6 ♦  10 

RETURN 

END 
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$  FASTRAN 

C  SUBROUTINE  DATEX3  6-8-67 

SUBROUTINE  DATEX3 

COMMON  BSCN, BSCN2 , BSFN , BSFN2  *  CPH *  CPHO , CTHO , DEN , EB , EBbC , EBSC2  * 

1  EM IN,ESO,EXC, FNDAX, IHIX, IHP^, I  RA  . I RC « LOOP , MXAX *  NDAX ♦ NE ,  NESO  , 

2  NESX«NMAX,NMIN,NPAX,NRAN,NSCI  «  NXSO  «  NXSX  *  PCS  *  PHO  «  RAN  ,  s  »  i>PH  , 

3  SPHO ♦ ST HO ,  SUR V  *  SUZ ,  TESC,  TESC2 ,  TFEX , TFEX2  ,  TF2N , TF2N2 » 1 SFN ♦ T SFN2 < 

4  TS 1 N« TS 1 N2  « TS2N ,  TS2N2 *  UT  ,  UX  ,  UY  »  UZ  * UZC T ♦ W  ,  WMAX • W I  , W I SO  * W I SOC  * 

5  X , XN , XR , XRS , XR2  *  XS , XS2 , XSECB , XSO  t Yt YN  <Zi ZN , ZP «  ZR ,  ZS 
DIMENSION  BSCN ( 30) ,BsCN2(3u) * BaFNi 3u ) *BoFN2t3u) »CPH,36o) , 

1  EB ( 25 ) «  EBSC ( 30 ) ,EBSC2(30) *EbO(30) ,EXC(899) , PCS (404) ,SPH(360> , 

2  SURV ( 899 ) »  TESC ( 30 i 30 ) »  IESC2(30,30) » lSFNt3o,3u) ,  lsFN2i3o»3u) , 

3  TS1N(30»30) »  TS  1 N2 (30*30) ,UT<899) ,WTS0(30) ,WTS0C(30) , 

4  XSECB (25,4) «  XSO  ( 80 )  »  ZP ( 30 ) 

DIMENSION  XSEC ( 4 ) 

NT  AB I N= 1 

DO  30  M= 1  «  MXAX 

SURV8= ( XSECB ( M « 1 )+XSECB(M,2) )/XSECB(M,3> 

PCSB=XSECB(M, 1 ) /(XSECB (M, 1 }+XSECB(M,2) ) 

XSECB ( M , 1 ) =SURVB 

XSECB ( Mi  2 ) =PCSB 

XSECB ( M,  3 ) =XSECB ( M , 3 ) *DEN 

DO  20  1=1*4 

20  XSECB (Mil ) =ELOG ( XSECB l  M ,  I  )  ) 

30  EB( M) =ELOG( EB ( M  )  ) 

MXAX2=MXAX/2 
MM=MXAX+ 1 
DC  40  M=1  ,  MX AX 2 
MM=MM— 1 
DUM=EB(M) 

EB(M)sEB (MM) 

EB(MM) =DUM 
DO  40  I =1 i4 
DUM= XSECB ( Mi  I ) 

XSECB (M« I  ) aXSECB ( MM  *  I ) 

40  XSECB (MM, I ) =DUM 
DO  60  NE= 1,899 
FNE=NE 

E=1000.0/( FNE+90 • 5 ) - 1 .0 
E=ELOG(E) 

CALL  T AB I N ( NT  AB I N , XSECB « EB , MXAX , 4 , E , XSEC ) 

SURV (NE)=EXP( XSEC ( 1 ) > 

UT  ( NE ) =EXP ( XSEC ( 3 )  ) 

EXC ( NE ) aEXP ( E+XSEC ( 4 ) ) 

I F ( NE— 404 )  50,50,60 
50  PCS ( NE )=EXP(XSEC(2> ) 

60  CONTINUE 
PH=-0  •  5 

DO  70  I PH= 1 ,360 
PH=PH+1 ,0 

PHR=PH  *0.017453293 
CPH ( 1PH) =COS(PHR) 

70  SPH(  I  PH ) =S I N (  PHR ) 

FNDAX=NDAX 
XR2=2.0*XR 
XS2=2 . 0*XS 
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WMAX=0.51 1 /EMIN 
FNXSX=NXSX 
DX=XR2/FNXSX 
XX=  — XR— DX/2  4 0 
DO  SO  NXSO=l « NXSX 
XX*XX+DX 
80  XSO ( NXSO ) a  XX 
I HPP= I H I X/NXSX 

I F ( NXSX* I HPP J I  H  I X )  90.120*90 
90  WOT  6.100 

100  FORMAT (61 H  NUMBER  OF  HISTORIES  NOT  DIVISIBLE  BY  NUMBER  OF  SOwRGE  P 
10INTS) 

110  CALL  SYSTEM 
120  IRC=NRAN 

STHO=SQRT ( 1 . 0-CTH0**2 ) 

PHO=PHO*0  .0 1 7453293 
CPHO«COS ( PHO ) 

SPHO=SIN(PHO> 

TS2N=0*0 
TS2N2=0.0 
TF2N-0.0 
TF2N2*0*0 
TFEX*0*0 
TFEX2=0.0 
DO  130  ND= 1 ♦ NDAX 
BSCN(ND) =0.0 
BSCN2(ND) =0.0 
BSFN(ND) =0  a  0 
BSFN2(ND1 =0.0 
EBSC( ND) *0.0 
EBSC2 ( ND ) =0  4  0 
DO  130  NP= 1 »  NPAX 
TS1N(NDoNP)*Q40 
TS1N2 (ND.NP) *040 
T C  "M  (  ND  .  NP  1  «*  0  i  0 
TSFN2(ND.NPI«0*0 
TESC ( ND . NP l«0.0 
130  TESC2(ND»NP)  a»0*  0 
WTSOC ( 1 ) =  WTSGH  1  ) 

IF(NESX“1 ) 160. 160. 140 
140  DO  150  NES0*2«NESX 

150  WTSOC (NESO)** WTSOC (NESO-1 )+WTS0<NE3Q) 

160  DO  170  NESOfe 1 t NESX 

WTSO ( NESO  J  *  WTSO  <  NE SO ) /WTSOC  C  NESX ) 

1  70  WTSOC ( NESO ) * WTSOC ( NESO ) /WTSOC  <  NESX  > 

RETURN 

END 
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FASTRAN 

SUBROUTINE  TRAGK3  6-7-67 

SUBROUTINE  TRAGK3 

COMMON  BSCN,  BSCN2  ,  BSFN  ,  BSFN2  ,  CPH*  CPHO ,  CTHO  *  DEN  4 EB  ♦  EBbC*  EB^C2  4 
EM  IN*  ESO  « EXC * FNDAX « IHIX* IHPP* IRA* I RC « LOOP *MXAX * NDAX *NE 4 NES6 « 


NESX «  NMAX  ♦ NM I N , NPAX *  NRAN « NSCT  «  NXSO * NXSX * PCS  * PHO  4  RAN * S  4  SPH  4 
SPHO*  STHO*  SURV  »  SUZ « TESC*  TESC2  *  TFEX  * TFEX2  *  TF2N*  TF2N2  «  TSFN«  JSFN2  4 
TSl N « TS1 N2  * T  S2N  «  TS2N2 « UT *  UX  *UY iUZi UZCT *W«WMAX«WT  « WTSQ*  WTSOG  « 


X « XN*  XR « XRS  ♦ XR2 1XS1 XS£  % XSECB  «  XSO  •  Y  «  Y!  1 »  Z  « ZN * ZP « ZR *  ZS 
DIMENSION  BSCN( 30)  *  BSCN2  <  30 ) *BSfN(30> *BSFN2<30 ) « CPH (360) 4 

EB( 25 ) « EBSC ( 30 ) « EBSC2( 30 ) *  ESO( 30) *EXC  ( 899) * PCS C  404 ) ♦ SPH ( 360 ) 4 
SURV ( 899 ) * TESC (30*30) * TESC2 ( 30  *  30 ) *TSFN(30*30) * TSFN2 ( 30 4 30 ) 4 
TSl N(30*30 ) *  TS1N2 (30*30 ) *  UT  <899 )  * WTS0(30) ♦ WTSOG (30 ) 4 
XSECB ( 25  *4 ) ♦ XSO (80)  *ZP(30) 

CALL  RANDCt IRC 4  RAN) 

IRA* IRC 
X=XSO(NXSO) 


Z  =  0.0 


Y=0«0 


SUZ=STHO 

UZ*CTHO 

UX=STHO*CPHO 

UY=STHO*SPHO 

CALL  RANDA (  I RA  <  RAN ) 

DO  4  NES0=1 *  NESX 
I F ( WTSOC ( NESO) —RAN ) 4  *  4  *2 
2  W=0 . 5 1 1 /ESO ( NESO ) 

GO  TO  6 
4  CONTINUE 
6  WT=1.Q 

NE= 1 000  « 0/ ( (0.51 1/W)  +  1 .0  )-90*0 
N=0 

CALL  RANDA (  I RA  « RAN ) 

S=— ELOG(  RAN ) /UT ( NE ) 


XN*X+S*UX 
YN*Y+S*UY 
ZN=Z+S*UZ 
CALL  CHECKS 


X=XN 

Y=YN 

Z*ZN 

IF(NSCT)  70* 10*2° 

N=N+1 

I  F  (  N—NMAX )  30*30*90 
NE=NE 

WT* WT*SURV ( NE ) 

IF(NE-404)  40*40*60 
CALL  RANDA  (  1  RA  ♦  RAN  ) 

I F ( RAN— PCS ( NE ) )  60*50*50 
WT=WT*2.0 

CALL  RANDA (IRA*  RAN ) 

UZ=-1 »0+RAN*2  *0 
SUZ=SQRT I  1 t0-UZ**2> 

CALL  RANDA < IRA  *  RAN) 


NPH=36C«0*RAN 
UX=SUZ*CPH(NPH+1  ) 
UY=SUZ*SPH < NPH+  1  ) 
W=1  .0 
NE=571 
GO  TO  10 
CALL  C0MPT3 
IF<W-WMAX)  1 0,90»90 
IF(N-NMIN)  90 ♦ 80  *  80 
CALL  GRADE3 
RETURN 
END 


(ft  u 
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FASTRAN 

SUBROUTINE  CHECK3  7-18-67 

SUBROUTINE  CHECKS 

COMMON  BSCN. BSCN2 , BSFN , BSFN2 ,  CPH ,  CPHO , CTHO  , DEN , EB  * EBSC . EBSC2  * 

1  emin.eso.exc.fndax, ihix, ihpp, ira. irc.loop.mxax.ndax.ne.neso. 

2  NESX . NMAX , NM I N . NPAX , NRAN. NSCT ♦ NXSO. NXSX . PCS , PHO ♦ RAN , S . SPH  , 

3  SPHO . STHO, SURV , SUZ *  TESC.  TESC2 , TFEX, TFEX2 *  TF2N. TF2N2, TSFN. TSFN2 ♦ 

4  TS 1  N » TS 1  N2 1  TS2N  *  TS2N2 ,  UT  .  UX  .  UY  *  UZ  *  UZCT  *  W  ,  WMAX  ♦  WT  «  W  I  SO ,  W  I  S.OC « 

5  X ♦ XN, XR ♦ XRS , XR2 » XS i XS2 *  XSECB  *  XSO  *  Y » YN  t Z  *  ZN «  ZP ♦ ZR . ZS 
DIMENSION  BSCN( 30) »BSCN2(30) .BSFN(30) «BSFN2(30 ) ♦ CPH (360) « 

1  EB ( 25 ) ♦ EBSC <  30) «EBSC2(30) *ES0(30) iEXC(899) » PCS (404) .SPH(360) » 

2  SURV  ( 899  >  ♦  TESC  ( 30  *  30  )  .  TESC2  (  30  *  30  )  ♦  TSFSM  (  30 . 30  )  ♦  TSFN2  (30.30)* 

3  TS1N( 30.30) . TS1N2( 30*30 ) .UT(899 ) .WTS0(30) . WTS0C(30) * 

4  XSECB ( 25  *  4 ) . XSO ( 80 )  . ZP ( 30 ) 

1F(UX) 10.290.20 

10  XN=  — XN 
PUX=-UX 
GO  TO  30 
20  PUX=UX 

30  I F ( ZN— ZR ) 40 . 280  «  280 
40  I F ( ZN ) 230 .230.50 
50  IF(ZN-ZS)2lO, 170,60 
60  IF(XN+XS) 140.70 .70 
70  1 F ( XN-XS ) 80 . 100 . 130 
80  ZN=ZN+UZ*(XS-XN)/PUX 
XN=XS 

I F ( ZN— ZR ) 90 « 120 .120 
90  IF(ZN-ZS) 1 10. 100. 100 
100  NSCT  =  0 

GO  TO  360 

110  XN=XN— PUX*(ZN-ZS)/UZ 
ZN=ZS 
GO  TO  100 

120  XN=XN-PUX*(ZN-ZR)/UZ 
ZN=ZR 
NSCT  =  —  1 
GO  TO  360 

130  IF(XN-XR) 140. 150. 160 
140  NSCT  = 1 

GO  TO  360 
150  XN=-XR 

GO  TO  140 
160  XN=XN— XR2 

I F ( XN-XR ) 60 . 150.160 
170  IF (XN-XR) 190. 150. 180 
180  XN=XN— XR2 
GO  TO  170 

190  IF (XN-XS) 200. 100, 140 
200  IF(XN+XS ) 140, 100, 100 
210  IF(XN-XR) 140, 150.220 
220  XN=  XN-XR2 
GO  TO  210 

230  XN=XN-PUX*ZN/UZ 
ZN=0  *  0 

240  1 F ( XN-XR )270,260, 250 
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250  XN=XN-XR2 
GO  TO  240 
260  XN=0«99999*XN 
270  NSCT=~1 

CO  TO  360 

280  XN=XN-PUX*(ZN-*2R)/UZ 
ZN«ZR 
GO  TO  240 

290  IF( UZ ) 340  *300  «  300 
300  IF(ABSF(XN)-XS) 310.310.330 
310  IF(ZN“ZS) 140.320.320 
320  ZN=*ZR 
NSCTa-1 
GO  TO  380 

330  IF(ZN-ZR) 140.320.320 
340  IF ( ZN 1 350 .3504 140 
350  ZN=0«  0 
NSCTs-1 
GO  TO  380 

360  1F( UX  > 370 . 380  «  380 
370  XN*~XN 
380  RETURN 
END 
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FASTRAN 

SUBROUTINE  C0MPT3  6~7-*67 

SAMPLE  NEW  DIRECTION  AND  ENERGY  FROM  COMPTON  DISTRIBUTION 
SUBROUTINE  C0MPT3 

COMMON  BSCN,  BSCN2 , BSFN , BSFN2 , CPH , CPHO , CTHO , DEN * EB  « EBSG  *  EBSC2  » 

1  EMIN,ESO,EXC, FNDAX, IHIX* 1HPP* IRA* I RC ,LOOP * MXAX * NDAX ♦ NE « NESO « 

2  NESX , NMAX  *  NM I N , NPAX , NRAN*  NSCT  « NXSO » NXSX . PCS ♦ PHO ♦ RAN ♦ S ♦ SPH  * 

3  SPHO , STHO , SURV , SUZ  *  TESC  *  TESC2 , TFEX  *  TFEX2 , TF2N  *  TF2N2 ♦ TSFN , TSFN2 ♦ 

4  TS1N«TS1  N2  <  T  S2N «  TS2N2  *  UT  *  UX ♦ UY « UZ*  UZCT  *  W ♦ WMAX  «  WT  *  WTSQ «  WTS&G  « 

5  X. XN ♦ XR *  XRS . XR2 *  XS , XS2 ♦ XSECB  *  XSO  *  Y «  YN ♦ Z *  ZN, ZP  *  ZR . ZS 
DIMENSION  BSCN ( 30)  ,BSCN2<30>  *BSFN(30) *BSFN2(30) *CPH(360) * 

1  EB ( 25 ) *  EBSC ( 30 ) ,EBSC2(30) ,ESO(30) *EXC(899) *  PCS (404) *SPH(360) 4 

2  SURV (899)  , TESC (30 *30 ) * TESC2 ( 30  *  30 ) * TSFN (30 *30)  *  TSFN2 ( 30 ♦ 30 ) 4 

3  TS1N(3C«30)  *TS1N2(30*30)  *UT(899)  *WTS0<30)  «WTSG)G(30)  4 

4  XSECB(25.4) «XS0(80) *ZP(30 j 
10  CALL  RANDA ( I RA * RAN ) 

T=2»0/W 

IF(RAN^( 1 .0+T)/(9.0+T) ) 20 *20 *30 
20  CALL  RANDA (  I RA  «  RAN ) 

R=  1 .0+RAN*T 

CALL  RANDA ( IRA* RAN) 

I F ( RAN— 4  .  0* ( R— 1  .0)/(R**2) ) 40* 40* 10 
30  CALL  RANDA <  IRA  4  RAN) 

R=  ( 1.0+T)/(1«0+RAN*T> 

CALL  RANDA ( IRA t RAN) 

I F ( RAN— 0 • 5* ( ( W  — R*W+ 1 *0 ) **2+ 1 .0/R) ) 40*40* 10 
40  WN=W*R 

COM=  1  •  0+ W-*WN 
W=WN 

I F ( W—WMAX )  45*80*80 
45  SOM=SQRT( 1 ,0-C0M**2) 

NE= 1 000 • 0/ ( (0.51 1/W)  +  1 «  0 ) u90  *0 
CALL  RANDA (  IRA  1  RAN) 

I PH  =  360  «  0*RAN 

UZN=UZ*COM+SUZ*SOM*CPH( IPH+1  ) 

SUZN=SQRT  (  1 . 0-“UZN**2  ) 

A=SUZ*SUZN 

IF ( A“0. 000001 >50,50*60 
50  UXN=-CPH ( I PH+ 1 ) *SUZN 
UYNsSPH (  I PH+ 1 ) *  SUZN 
GO  TO  70 

60  CDPH« (COM“UZ*UZN)/A 

SDPH°SOM*SPH ( IPH+1 ) /SUZN 
UXNs  <  <UX*CDPH-UY*SDPH)*SUZN)/SUZ 
UYNa ( <UY*CDPH+UX*SDPH) *SUZN j/SUZ 
70  UX=UXN 
UY=UYN 
UZ=UZN 
SUZ=SUZN 
80  RETURN 
END 
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FASTRAN 

SUBROUTINE  GRADE3  6-7-67 

SUBROUTINE  GRADE3 

COMMON  BSCN, BSCN2 ♦ BSFN . BSFN2 ,  CPH, CPHO , CTHO , DEN , EB  , EBSC »  EBSC2 ♦ 

1  EMIN,ESO,EXC,FNDAX, IHIX. IHPP. IRA, I RC , LOOP , MXAX ♦ NDAX , NE , NESO , 

2  NESX «  NMAX , NM I N , NPAX , NRAN, NSCT , NXSO  ,  NXSX ,  PCS  , PHO  .  RAN , S , SPH ♦ 

3  SPHO » STHO.  SURV , SU2 .  TESC . TESC2  rTFEX  * TFEX2 , TF2N , TF2N2  ,  TSFN , TSFN2 , 

4  TSl N. TS 1 N2 . TS2N , TS2N2 1 UT , UX , UY , UZ  * UZCT , W . WMAX » WT , WTSO , WTSOC  » 

5  X,XN.XR,XRS,XR2,XS,XS2,XSECB»XS0, YiYNtZ.ZN, ZP.ZR.ZS 
DIMENSION"  BSCN( 30} ,BSCN2(30) ,BSFN(30) ,BSFN2(30) ,CPH(360) , 

1  EB( 25) .EBSC (30) «EBSC2(30) .ES0C30) .EXC(899) .PCS (404) .SPH (360) . 

2  SURV (899) .TESC( 30.30) . TESC2 ( 30 . 30 ) .TSFN (30 .30) . TSFN2 ( 30 . 30 ) . 

3  TSl N( 30. 30) .TSl N2 ( 30 . 30 ) .UT(899) .WTSO (30) .WTSOC (30) . 

4  XSECB( 25 .4 ) . XSO ( 80 ) ,ZP(30) 

IF(UZ) 10, 10,20 

10  ND=(X+XR)*FNDAX/XR2 
XFN=— WT/UZ 
EX=EXC ( NE  5  *XFN 
BSCN ( ND+ 1 ) =8SCN ( ND+ 1 ) +WT 
BSCN2 ( ND+ 1 ) =BSCN2 ( ND+ 1 ) +WT**2 
3SFN( ND+1 )=BSFN(ND+1 )+XFN 
BSFN2 ( ND+ 1 ) =BSF  N2 ( ND+ 1  )+XFN**2 
EBSC ( NDTl  ) =EBSC ( ND+ 1 ) +EX 
EBSC2 ( ND+ 1 ) =EBSC2 ( ND+ 1 )+EX**2 
GO  TO  140 
20  WT2=WT**2 
XFN.WT/UZ 
XFN2»XFN**2 
EX»EXCCNE)*XFN 
EX2=EX**2 
TS2N»TS2N+WT 
TS2N2=TS2N2+WT2 
TF2N«TF2N+XFN 
TF2N2=TF2N2+XFN2 
TFEXaTFEX+EX 
TFEX2=TFEX2+EX2 
IF ( UZ—UZCT ) 140. 140,25 
25  DO  130  NP= 1 « NPAX 
SC= ( ZP ( NP ) -ZR ) /UZ 
XT=X+UX*SC 

IF( ABSF( XT) -XR) 90,30,40 
30  ND*0 

GO  TO  100 
40  IF(UX)70, 50,50 
50  NRW=XT/XR2 
FNRW=NRW 
XT=XT— FNRW*XR2 
I F ( XT —XR ) 90 , 30 , 60 
60  XT=XT“XR2 
GO  TO  90 
70  NRW=-XT/XR2 
FNRW=NRW 
XT=XT+FNRW*XR2 
IF(XT+XR>80, 30. 90 
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90  ND=<XT+XR)*FNDAX/'XR2 
100  TS1 N( ND+ 1 «  NP ) =TS 1 N ( ND+ 1 *NP»+WT 

TS1N2(ND+1 ,  NP) =TS 1 N2 ( ND+ 1 iNP)+WT2 
TSFN( ND+l  i NP ) rTSFN ( ND+ 1 4  NP ) +XFN 
T SFN2  <  ND+ 1 s  NP ) =  T SFN2 ( ND+ 1 «NP)+XFN2 
TESC ( ND+ 1 « NP ) = TESC ( ND+ 1 i NP ) +EX 
TESC2 ( ND+ 1 « NP ) =  TESC2 ( ND+ 1 « NP  } +EX2 
130  CONTINUE 
140  RETURN 
END 


o  » 
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FASTRAN 

SUBROUTINE  TEACH3  7-26-67 

SUBROUTINE  TEACH3 

COMMON  BSCN. BSCN2 * BSFN .  BSFN2  *  CPH . CpHO * CTHO . DEN ♦ EB ♦ EBSC , EBSC2 * 

1  •  emin,eso,exc,fndax, IHIX, IHPP* IRA* IRC* LOOP, mxax.ndax*ne*neso, 

2  nesx,nmax*nmin,npax*nran,nsct*nxso«nxsx*pcs,pho*ran*s,sph« 

3  SPHO , STHO  *  SURV , SUZ , TESC , TESC2  *  TFE X . TFEX2  *  TF2N , TF2N2 , TSFN , TSFN2 ♦ 

4  TS 1 N*  TS 1 N2  »  TS2N , TS2N2 , UT . UX , UV , UZ  *  UZCT ♦ W » WMAX , WT  * WTSO , WTSQG , 

5  X , XN*  XR *  XRS , XR2  *  XS  *  XS2  *  XSEGB *  XSO , Y , YN , Z  *  ZN , ZP  » ZR , ZS 
DIMENSION  BSCN( 30 ) *BSCN2(30) »BSFN(30) *  BSFN2  <  30 ) *GPH(36Q) * 

1  EB  <  25 ) , EBSC ( 30) *EBSC2<30) *ES0(30)  *EXC(899) *  PCS (404) ,SPH(36o>  « 

2  SURV (899) *  TESC ( 30  *  30 ) *  TESC2 ( 30 , 30  > , TSFN (30 *30) *  TSFN2 ( 30  *  30 ) * 

3  TSlN(30*30)  ,TSl N2 ( 30 ♦ 30 ) , UT ( 899 ) »WTS0(30) *  WTSQG  4  30 ) 4 

4  XSECB (25,4),XS0(80)  *  ZP ( 30 ) 

DIMENSION  TS3N( 30 ) *TS3N2<30) *TF3N(30) ,TF3N2(30) »TFEX3(30) ♦ 

1  T FEX32 ( 30 ) *  XDET ( 30 ) 

H I X= I  H  I  X 

FLOOPaLOOP 

H I X*H I X*FLOOP 

DO  40  NP= 1  *  NPAX 

TS3N(NP)a0*0 

TS3N2 ( NP ) *0 • 0 

TF3N ( NP |a0«0 

TF3N2 ( NP ) aO  «  0 

TFEX3 ( NP ) =0  4 0 

TFEX32 ( NP ) aO  *0 

DO  10  ND= 1 , NDAX 

TS3N <  NP )  =  TS3N ( NP ) +TS 1 N ( ND *  NP ) 

TS3N2 ( NP ) =TS3N2 ( NP ) +TS 1 N2  <  ND  *  NP ) 

TF3N( NP ) =TF3N ( NP ) +TSFN ( ND*  NP ) 

TF3N2 ( NP ) =  TF3N2  < NP ) +TSFN2  <  ND  *  NP ) 

TFEX3 ( NP ) =TFEX3 ( NP ) +TESC ( ND ♦ NP ) 

1 0  TFEX32 ( NP )  =  TFEX  32 ( NP ) +TESC2 ( ND  *  NP ) 

TS3N ( NP ) =TS3N ( NP  > /H I  X 
TS3N2(NP)=TS3N2 (NP)/HIX 
TF3N ( NP ) 2 TF3N ( NP ) /H I  X 
TF3N2  (  NP  )  =TF3N2  (  NP  )  /H I,X 
TFEX3 ( NP ) =TFEX3 ( NP ) /H IX 
TFEX32 ( NP )  a  TFEX32 ( NP ) /H I  X 
IF ( T  S3N( NP ) >30*20*30 
20  TS3N2(NP) a-lOiQ 
TF3N2  (  NP  ) =- 1 0  40 
TFEX32 ( NP ) au 1 0  4  0 
GO  TO  40 

30  TS3N2 ( NP ) =  SORT (  ( TS3N2 ( NP ) -TS3N ( NP ) **2 ) / ( H I  X- 1 . 0 ) )/TS3N(NP) 
TF3N2(NP)  =SQRT(  (TF3N2(NP)-TF3N(NP)**2)/'(HIX-1 .0)  >/TF3N(NP) 

TFEX32 ( NP \  =SQRT ( ( TFEX32 ( NP ) -TFEX3 ( NP ) **2 ) / ( H I X- 1 .0) )/TFEX3(NP) 

40  CONTINUE 

T S2N  =  TS2N/H I  X 
TS2N2  =  TS2N2/H  I  X 
TF2N=TF2N/HI X 
TF2N2=TF2N2/H I  X 
TFEX»TFEX/HIX 
TFEX2  =  TFEX2/H I  X 
I F ( TS2N) 60*50*60 
50  TS2N2=- 10*0 
TF2N2=— 1 0*0 


TFEX2=- 1 0  #  0 
GO  TO  70 

60  TS2N2=SQRT(  ( TS2N2-TS2N**2 ) / < H I X- 1 .0) ) /TS2N 
TF2N2=SQRT( < TF2N2-TF2N**2 ) / ( H I X- 1 . 0 > J/TF2N 
TFEX2=SQRT< ( TFEX2— TFEX**2 ) / ( HI XW1 .0) )/TFEX 
70  DO  130  ND= 1 4  NDAX 

BSCN<  ND) «BSCN( ND ) /H I X 
BSCN2 ( ND ) =BSCN2 ( ND ) /H I  X 
BSFN ( ND > =BSFN <  ND ) /H 1 X 
BSFN2 <  ND ) =BSFN2  <  ND ) /H I  X 
EBSC ( ND ) *EBSC ( ND ) /H I X 
EBSC2 ( ND ) =EBSC  2  <  ND ) /H I  X 
I F ( BSCN ( ND ) ) 90  «  80  «  90 
80  BSCN2(ND)=-10»0 
BSFN2 ( ND ) =— 10*0 
EBSC2 ( ND ) =— 1 0  #0 
GO  TO  100 

90  BSCN2 ( NO  >  =SQRT (  ( BSCN2 ( No ) -BSCN ( ND ) **2 > / 1 H I  X- 1 . 0 ) > KBSCN ( ND ) 
BSFN2  (  ND )  =SQRT  (  <  BSFN2  (  ND  )  -BSFN (  ND )  **2 )  / ( H I  X- 1  .0)  )  J^BSFNii ND  j 
EBSC2(ND)  =SQRT(  (  EBSC2  {  ND  )— EBSC  (  ND )  **2  )  /  (  HI  X“i  •  0)  )V%BSc'<NDi 
100  DO  130  NP=1«NPAX  ' 

TS 1  N  (  ND  *  NP  )  =:  TS 1  N  (  ND  «  NP  )  /H I X 
TS1N2 (NDtNP) =TS1N2<ND,NP)/HIX 
TSFN ( ND »  NP ) = TSFN ( NDiNP)/H1X 
TSFN2 <  ND . NP )  =  T  SFN2 ( ND , Np ) /H I  X 
TESC (ND»NP)«TESC(ND.NP)/HIX 
TESC2<ND«NP) =TESC2(ND«NP)/hlX 
IF(TS1N(ND,NP) > 120* 1 10 » 120 
110  TS1N2<ND*NP)=~10.0 
TSFN2 ( ND  «  NP ) =-l0.0 
TESC2  (  ND  •>  NP )  =«*  1  0  «0 
GO  TO  130 

120  TS1N2(ND.NP)=SQRT< < TS 1 N2 < ND » NP ) -TS 1 N < ND» NP ) **2 > / < H1H- 1 .0) ) 
1  /TS1 N ( ND  *  NP )  1  l 

TSFN2(ND,NP)=SQRT( (TSFN2{ND.NP)-TSFN<NDfNP)**2)/(HlX“l .0) ) 

1  /TSFN(ND4NP)  ’  '  I  ' 

TESC2 ( ND  s  NP ) =SQRT(  ( TESC2 < ND . Np ) -TESC < ND . NP ) **2 ) ✓ ( H I X- 1 • 0 ) ) 

1  /TESC(ND4NP)  '  I 

130  CONTINUE 
D0SO=»0*0 

DO  140  NES0=1 4 NESX 
NE= 1 000 «  0/ <ESO ( NESO )  + 1 «0)“90i0 
140  DOSO^DOSO+EXC  <  NE ) *WTSO  <  NESO ) 

TF2N«TF2N*CTHQ 
TFEX=TFEX*CTHO/DOSO 
DO  150  NP= 1 i NPAX 
TF3N <  NP ) =TF3N ( NP  >  *CTHO 
1  50  TFEX3 ( NP ) =TFEX3 ( NP ) *CTHO/DOSO 
DO  160  ND=»  1  i  NDAX 
BSCN  (  ND )  siBSCN  (  ND  >  *FNDAX 
BSFN<  ND) »BSFN<  ND)*FNDAX*CTHO 
EBSC  (  ND )  tiEBSC  <  ND  )  *FNDAX*GTHO/DOSQ 
DO  160  NPnl , NPAX 
TS 1 N ( ND  *  NP ) =TS 1 N ( ND  « NP ) *FNDAX 
TSFN ( ND . NP ) d TSF N  <  ND  »  NP ) *FNDAX*CTHO 
1  60  TESC  <  ND » NP )  =  TESC ( ND . NP ) *FNDAX*GTHO/DOSO 
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DO  16.5  NP=  1  «  NPAX 

165  2P(NP).'sZP<  NPl^ZR 

170  FORMAT (1 HO) 

175  FORMAT ( 1H1 ) 

W0T6«175  f 

WOT 6  4 1 70 
W0T64l80< IRC 

180  FORMAT ( 21 H  FINAL  RANDOM  NUMBER=013) 

WOT  6 4 170 
WOT 6 ♦ 1 85  4DOSO 

185  FORMAT (54H  EXPOSURE  CURRENT  DUE  TO  UNIT  INCIDENT  NUMBER  CURRENT = 

1  1PE10.2) 

W  OT  6  4  1 7  0 

DX=XR2/FNDAX 

XX=-DX/2.0 

DC  190  ND=  1  4  NDAX 

XX=XX+DX 

190  XDET ( ND ) =XX 
W0T6«200 

200  FORMAT (32H  BACKSCATTER I NG  FROM  RIBBED  SLAB) 

WOJ6«210 

210  FORM AT ( 72H  NUMBER  CURRENT  (FLUX)  NORMALIZED  TO  UNIT  INCIDENT  NUMBE 
1R  CURRENT  (FLUX)) 

W0T6 1 220 

220  FORMAT ( 46H  EXPOSURE  NORMALIZED  TO  UNIT  INCIDENT  EXPOSURE) 

W0T6. 170 
WOT  6 1 230 

230  FORMAT (11H  HORIZONTAL) 

WOT  6 4 240 

240  FORMAT ( 31 H  DETECTOR  NUMBER  NUMBER) 

W0T6 4 250 

250  FORMAT ( 4 1 H  POSITION  CURREN t  FLUX  EXPOSURE) 

DO  260  ND=1 4  NDAX 

260  W0T6  «270  4  XDET ( ND ) *  B5CN( ND ) 4 BSFN ( ND > 4 EBSC ( ND ) 

270  FORMAT (FI 103« 1P1 OE 10.2) 

W0T6  4 175 
W0T6« 170 
WOT 6 « 280 

280  FORM AT ( 64 H  FRACTIONAL  STATISTICAL  DEVIATION  OF  RESULTS  IN  PRECEEDI 
1 NG  TABLE) 

WOT6o 170 
WOT 6 « 230 
WOT  6 » 240 
WOT  6  4  250 

DO  290  NO= 1 9  NDAX 

290  WOT6  4300  4  XDET ( ND ) 4  BSCN2 ( ND ) 4  BSFN2 ( ND ) 4  EBSC2 ( ND ) 

300  F0RMAT(F1 1.34 10F10.4) 

WOT  6  « 175 
WOT  6  4 1 70 
WOT64310 

310  FORMAT (56H  TRANSMISSION  AVERAGED  OVER  ALL  DETECTORS  IN  GIVEN  PLANE 
1  ) 

WOT 6  4  320 

320  FORMAT ( 1 7H  NO  COSINE  CUTOFF) 

W0T64330 

330  FORMAT ( 48H  RESULTS  NORMALIZED  AS  IN  CASE  OF  BACKSCATTER I NG ) 
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WOT  6  < 170 
WOT 6 *340 

340  FORMAT (371'  FRACT 

WOTfc . 350 

350  FORMAT  ( 3? K  .-T  L^V.WT 

W0T6*360*TS2N*  TS2N2 

360  FORMAT  (  1 6H  NUMBER  CURRENT  1PE1 0. 2 *  3,"'F1 0. 4 ) 
W0T6  ♦ 370 , TF2N*  TF2N2 

370  FORMAT ( 1 6H  NUMBER  FLUX  1  PE  1 0 . 2 , OPFX 0 . 4 ) 
W0T6 « 380 «  TFEX  *  TFEX2 

380  FORMAT ( 1 6H  EXPOSURE  1  PE  1 0 . 2  * 0PF1 0 . 4 > 

W0T6* 175 
W0T6* 170 
W0T6*310 
WOT 6  *  390  « UZCT 

390  FORMAT ( 1 5H  COSINE  CJT0FF=F965) 

WOT 6 * 330 
W0T6  *  1 70 
W0T6  *  400 

400  FORMAT ( 1 1 H  HEIGHT  OF) 

WOT  6*410 

410  FORMAT ( 1 1 H  DETECTOR) 

W0T6*420 

420  FORMAT ( 31 H  PLANE  NUMBER  NUMBER) 

W0T6  *  430 

430  FORMAT ( 4 1 H  ABOVE  RIB  CURRENT  FLUX 


FRACTIONAL) 


HEIGHT  OF) 


DETECTOR) 


PLANE 


NUMBER 


430  FORMAT ( 4 1 H  ABOVE  RIB  CURRENT  FLUX  EXPOSURE) 

DO  440  NP=1*NPAX 

440  W0T6*270*ZPCNP) ,TS3N(NP) ,TF3NCNP) ,TFEX3(NP) 

WOT  6 ♦ 1 75 
WOT  6  *  1 70 
WOT  6 « 280 

WOT  6*1 70  „ 

WOT  6 « 400 
WOT  6*410 
WOT 6  *  420 
WOT  6  *  430 

DO  450  NP=  1  *  NPAX 

450  WOT 6*300*ZP( NP )  . TS3N2 ( Np ) , TF3N2 < NP ) *TFEX32<NP) 

W0T6* 175 
WOT  6  *  170 
WOT 6 « 460 

460  FORMAT! 100H  IN  THE  TABLES  WHICH  FOLLOW*  H=  HEIGHT  OF  DETECTOR  PLAN 
IE  ABOVE  RIB*  X=  POSITION  OF  DETECTOR  IN  PLANE) 

W0T6* 170 
WOT 6  *470  *  UZCT 

470  FORMAT ( 43H  TRANSMITTED  NUMBER  CURRENT*  COSINE  CUT0FF=F9.5) 

WOT 6  *  480 

480  FORMAT (43H  NORMALIZED  TO  UNIT  INCIDENT  NUMBER  CURRENT) 

CALL  R I BTAB  <  NPAX  *  NDAX *  ZP • XDET , TS 1 N ) 

WOT 6  *  175 
W0T6* 170 
WOT  6 ♦ 280 

CALL  R I BT AB ( NPA X  *  NDAX  *  ZP ♦ XDET ♦ TS 1 N2 ) 

WOT  6  *  175 
WOT 6  *  170 
WOT  6*490*  UZCT 


NUMBER) 


FLUX  EXPOSURE) 


490  FORMAT ( 40H  TRANSMITTED  NUMBER  FLUX*  COSINE  CUTQFF=F9.5) 

WOT 6 « 500 

500  FORMAT C40H  NORMALIZED  TO  UNIT  INCIDENT  NUMBER  FLUX) 

CALL  R I  STAB ( NPAX  . NDAX *  ZP ♦ XDET , TSFN ) 

W0T6< 175 
W0T6. 170 
W0T6.280 

CALL  RIBTAB(  NPAX. NDAX. ZP.XDE1'.  I  SFN2  ) 

W0T6« 175 
W0T6« 170 
W0T6.510.UZCT 

510  FORMAT (64H  EXPOSURE  DUE  TO  iRANbMIllED  bCA I  I ERED  RAD  I A i  I  ON.  COblNE 
1  CUTQFF=F9*5) 

W0T6«520 

520  FORMAT (54H  NORMALIZED  10  UN  I  I  EXPOSURE  DUE  |0  INCIDENi  RAD I A  I  1  ON ) 
CALL  RIBTABINPAX.NDAX.ZP.XDE  I  . t ESC  )  1 

WQT6.175 
W0T6. 170 
W0T6.280 

CALL  RIBTAB(NPAX.NDAX«ZP<XDET«TESC£) 

WOT  6 « 175 

W0T5.530. (ZP(NP) ,NP=1 .NPAX) 

530  FORMAT! 1P7E1 0.3) 

W0T5.530, (XDET ( ND) ,ND=1 .NDAX) 

W0T5.530. ( ( TESC  (ND.NP) ,ND=1 .NDAX) .NP=1 .NPAX) 

WOT 5 . 530 . ( (TESC2( ND.NP) ,ND=1 .NDAX) «NP=1 ,NPAX) 

NT  AB I N= 1 

DO  540  NPa 1 . NPAX 
ZT*ZP ( NP ) +ZR 
DO  540  ND= 1 . NDAX 

TESC2  (ND.NP)  =>TESC2  ( ND  .  NP  )  *TESC  (ND.NP) 

CALL  DIRIB(NTABIN»EB,XSECB( 1 .3) . MXAX. CPHO. C I  HO ♦ XR . Xb , /R . ZS ♦ 

1  XDET(ND).ZT  .ESO.WTSO.NESX,TSFN(ND<NP) ) 

TSFN  (  ND  «  NP  )(*  TSFN  (  ND  .  NP  )  *eTHQ 

TESC (ND.NP) a  TESC (ND.NP) +TSFN (ND.NP) 

540  TESC2(ND.NP) =TESC2( ND.NP) /TESC (ND.NP) 

W0T6. 170 
WOT  6 « 550 

550  FORMAT ( 38H  EXPOSURE  DUE  TO  UNSCATTERED  RADIATION) 

WOT 6 . 520  1 

CALL  R I BTAB  (  NPAX .  NDAX .  ZP .  XDET  ,  JTSFN  ) 

WOT  6.1 75 
W0T6. 170 
WOT 6 . 560 . UZCT 

560  FORMAT (54H  EXPOSURE  v  SCATTERED  PLUS  UNSCATTERED > .  COSINE  CUTOFF* 

1  F9.5I  I  'I 

WOTS .520 

CALL  R IBTAB (NPAX. NDAX. ZP. XDET. TESC ) 

WOTS .175 
WOT 6.170  * 

WOT  6 . 280 

CALL  R I BT AB ( NPAX . NDAX , ZP » XDET . TESC2 ) 

WOT  6 . 175 

W0T5.530. ( (TESC  (ND.NP) .ND=1 .NDAX) *NP=1 .NPAx) 

RETURN 


n  w 
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fastran 

SUBROUTINE  DIR  IB  7-1-67 

SUBROUT  I NE  D I R I B ( NTA6 IN*  EB , XSECB , MXAX  , CPHO  , CTHO , XR , XS , ZR , ZS , XD , ZP  * 
1  ESO,WTSO*NESO,DIR) 

DIMENSION  ESC  25)  , XSECB (25,2) ,ESO(30) ,WTS0<30) , UAT (30) ,UEN(30) ♦ 

1  XSEC ( 2 ) 

UX= CPHO* SORT { 1 . 0-CTH0**2 ) 

IF(UX) 1 0 i 20 <  20 
10  UX=-UX 
X=XR— XD 
GO  TO  30 
20  X=XD— XR 
30  T  =  0  » 0 
ZT=ZR 

X=X-UX* ( ZP-ZR ) /CTHO 
aO  IF(  ABSF(X) -*XR)  70*60  4  50 
50  X=X+2«. 0*XR 
GO  TO  40 
60  X=XR 

70  IF(ABSF(X)**XS)  130*80*80 
80  IF(X“XS)90, 130, 100 
90  X=X+2«0*XR 
100  XT=X-UX*<ZT-ZS) /CTHO 
IF(XT-XS) 120,110,110 
110  T=T+ZT/CTHQ 
GO  TO  160 
120  S=<X-XS)/UX 
T=T+S 

ZT=ZT-CTHO*S 

X=XS 

130  XT=X-UX*(ZT-ZS)XCTHO 
I F ( XT+XS )  150, 140, 140 
140  T=T+ZS/CTHO 
GO  TO  160 

150  ZT=ZT~CThO*(X+XS)/UX 
X=2,0*XR-4XS 
GO  TO  100 

160  GO  TO  < 170, 190)  4NTABIN 
170  DNORMff 0  »Q 

DO  180  NES=1  4  NS SO 
EL^ELOGCESOINES) ) 

CALL.  TAB  I  N  { NT AB  I N ,  XSECB  ,  EB ,  MXAX ,  2  4  EL  *  XSEC  ) 

UAT ( NES ) =EXP ( XSEC ( 1 )  ) 

UEN ( NES ) =EXP I XSEC ( 2  >  +EL  ) 

180  DNORM^DNORM+WTSO ( NES ) *UEN  <  NES ) 

190  D I R*0  *  0 

DO  200  NES= 1 , NESO 

2  00  D I R=D I R  + WTSO ( NES  >  *U£M ( NES ) *EXP ( -UAT ( NES ) *T ) 

DIR=»DIR/<DNORM*CTHO) 

RETURN 

END 
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SUBROUTINE  RIBTAB  5-22-67 

SUBROUT  I NE  R I BT AB ( NZAX » NXAX  *Z»X» A) 
DIMENSION  ZOO)  »  X  (  30  )  ♦  A  (  30 »  30  ) 

W0T6  ♦  1  0 
10  FORMAT! 1 HO) 

I F ( NZAX- 1 0 ) 30  ♦  20  « 20 
20  NZT  =10 
GO  TO  40 
30  NZT =NZAX 
40  NZB= 1 

50  W0T6  <60 ♦ ( Z ( NZ ) «  NZ=N2B  *  NZT ) 

60  FORMAT ( 1 2H  X*  H=F7< 3 « 9F 1 0 . 3 ) 

DO  70  NX= 1 «  NXAX 

70  W0T6«80»X( NX)  «  < A(NX«NZ) «NZ=NZB*NZI  ) 

80  FORMAT (  F 1 0  •  3  » 1P1 OE 10*2) 

W0T6(- 10 

IF ( NZT— NZAX) 90  « 1 10- '  10 
90  NZB=NZT+1 
NZT=NZT+10 

I F ( NZT— NZAX ) 50  «  50  « 1 00 
100  NZT=NZAX 
GO  TO  50 
1 10  RETURN 
END 


O  tfl 
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FASTRAN 

SUBROUTINE  TAB  IN  28-8-64 

SUBROUT  I NE  T AB I N ( NTAB I N , FB . XB , NMAX » MM AX ♦  X ♦ FX 1 
DIMENSION  FB<  25*8) , XB ( 25 ) iFX ( 8 ) ,XBAV(25) »D 1( ^5 ) , D2 < 25 ) 
GO  TO  ( 10,30) iNTABIN 
10  NMAX1 =NMAX— 1 

DO  20  N=2,NMAX1 

X8AV { N ) = ( XB ( N- 1 ) +XB ( N ) ) /2 . 0 

D 1  (  N  )  =  ( XB ( N— 1 ) — XB ( N ) )*CXBCN-1 )-X6CN+l )  ) 

20  D2 ( N )  =  ( XB ( N ) —X  B  <  N— 1 ) ) * ( XB ( N ) -XB ( N+ 1 ) ) 

NT  AB I N=2 
NEXS=NMAX 1 —2 
30  IF(X-XB( 1 1 )60,50,40 
40  NX=2 

GO  TO  200 
50  NX=1 

GO  TO  220 

60  IF(X-XB(2) 180,70,40 
70  NX=2 

GO  TO  220 

80  IF(NEXS)90, 1 10, 140 
90  W0T6 , 1 00 

100  FORMAT (33H0N0T  ENOUGH  BASE  POINTS  FOR  TAB  I N ; 

CALL  SYSTEM 

110,  IF(X-XB(  NMAX  )  )  130,  120, 130 
120  NX=NMAX 

GO  TO  220 
130  NX=NM AX  1 
GO  TO  200 

140  DO  170  N=3 , NMAX 1 

I F ( X-XH ( N ) ) 1 70 , 1 50 , 1 60 
150  NX-N 

GO  TO  220 
160  NX=N 

GO  TO  180 
170  CONTINUE 
'  GO  TO  110 

180  IF(X-XBAV(NX> 1200,200, 190 
190  NX=NX— 1 

200  WT 1 = ( X-XBC  NX) ) * (X-XB(NX+1 )  ) /D l (NX) 

WT2=  <  X-XB ( NX-1 )  1  * ( X-XB ( NX+1 1 1 /D2 ( NX ) 

WT3= 1 • O-WT 1  — WT  2 
DO  210  M= 1 , MMAX 

210  FX(M)=WT1*FB(NX-1 ,M)+WT2*F3(NX,M)+WT3*FB(NX+1 ,M) 

RETURN 

220  DO  230  M= 1 , MMAX 
230  FX(M) =FB (NX, M) 

RETURN 

END 
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RANDA 

RANDB 

RANDC 

RANOD 


MULT 

TEMP 


SCATRE 

.PUNCH  OBJECT 

REM 

CALL  RAND  A.B.C.D 

< IR.R) 

REM 

FIRST  DIGIT  OF  IR 

LESS  OR  EQUAL 

3  OCTAL 

REM 

LAST  DIGIT  OF  IR  EQUAL  1  BINARY 

REM 

AT  LEAST  7  DIGITS 

TO  EXPRESS  IR 

IN  OCTAL 

ENTRY 

RANDA 

ENTRY 

RANDB 

entry 

RANDG 

ENTRY 

RANDD 

CL  A 

MULT 

TRA 

*+6 

cla 

MULT+1 

TRA 

*+4 

CLA 

MULT +2 

TRA 

*+2 

CLA 

MULT+3 

STO 

Temp 

LDQ* 

1  44 

MPY 

TEMP 

stq* 

1  44 

PXD 

0,0 

LLS 

27 

ADD 

*+4 

FAD 

*+3 

STO* 

244 

TRA 

344 

OCT 

200000000000 

OCT 

343277244615 

OCT 

011060471625 

OCT 

000272207335 

OCT 

000007346545 

BSS 

1 

END 
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